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Abstract 

a; 

In an earlier publication, we have introduced a method to obtain, at large N, the effective action for 
5/3 ' d-dimensional manifolds in a iV-dimensional disordered environment. This allowed to obtain the Func- 

tional Renormalization Group (FRG) equation for JV = oo and was shown to reproduce, with no need 
for ultrametric replica symmetry breaking, the predictions of the Mezard-Parisi solution. Here we com- 
pute the corrections at order 1/N. We introduce two novel complementary methods, a diagrammatic 
and an algebraic one, to perform the complicated resummation of an infinite number of loops, and derive 
the /3-function of the theory to order 1/N. We present both the effective action and the corresponding 
functional renormalization group equations. The aim is to explain the conceptual basis and give a de- 
tailed account of the novel aspects of such calculations. The analysis of the FRG flow, comparison with 
other studies, and applications, e.g. to the strong-coupling phase of the Kardar-Parisi-Zhang equation 
i—i ' are examined in a subsequent publication. 

>' 

»; 

^ : 1 Introduction 

In a series of recent articles we have constructed the Functional Renormalization Group (FRG) method 
for disordered systems, applied to specific situations and beyond one loop [1-9]. This method is, apart 
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from mean field theory [11,12] using Replica Symmetry Breaking (RSB) and some rare exactly solvable 
cases [13,14], the only known analytical method which promises to handle the strong coupling glass phase 
of disordered elastic systems [1-9,15-38]. Such systems, modeled by an elastic manifold (of internal 
dimension d) with a iV-component displacement field u(x) (i.e. x G M d and u(x) G W N ), are of high 
interest for numerous experiments [39-43,25,36,44]. This so-called random manifold model still offers 
great theoretical challenges and a strong motivation is the hope to gain insight into glassy physics. In 
addition, the d = 1 case maps onto the much studied Kardar-Parisi-Zhang growth equation [45]. It exhibits 
^ a strong coupling phase for which the upper critical dimension is still under debate [46-50]. 

Higher loop studies of the statics of disordered elastic systems allow, in principle, a systematic di- 
mensional expansion, in the simplest case around d = 4. They are however of a rather different nature 
than in standard field theory for pure critical systems [1,7,10,31,33-35,51-54]. Thermal fluctuations are 
found to be formally irrelevant in these glass phases, suggesting that the physics is controlled by a zero 
temperature fixed point. However before this fixed point is reached, the zero temperature effective action 
is found to become non-analytic [17]. Although this allows to evade the so-called dimensional reduction 
[55] which makes naive perturbation theory useless and yields unphysical results, it also generates amazing 
new subtleties in the field theory. These were analyzed in a number of papers [1,7,10,31,33-35,51-54], 
and although some solutions to the puzzles were proposed the physics still remains to be elucidated. 
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An interesting limit where one can hope to gain insight into these formidable problems is the large- 
TV limit. Since iV = oo is formally the mean-field limit, it allows a direct confrontation between the 
FRG method and mean field methods. A solution of the O(N) random manifold model for N = oo was 
proposed by Mezard and Parisi, using a saddle point with spontaneous replica symmetry breaking [11,12]. 
As in other models of glasses, spontaneous RSB can be related to ergodicity breaking of the Gibbs measure 
into several ground states [56]. Although it offers a rather elegant way out of dimensional reduction, it is 
by no means clear that systems with (large but) finite N should exhibit such a tremendous degeneracy of 
low energy states; and there are in fact indications to the contrary [57]. 

It is thus crucial to develop another line of attack, even in that limit. This is what we have achieved 
in a previous publication, where we have computed the effective action of the theory at large N. There 
we have derived the /3-function of the field theory to dominant order, i.e. for N = oo [2,3,8]. We have 
discovered that beyond the Larkin length the FRG flow freezes (at least for specific initial conditions) and 
that most of the features of the Mezard Parisi solution can be recovered. Interestingly however, in this 
formulation there is no need for a spontaneous RSB ansatz with ultrametric structure. Thus one may hope 
that it could be more adapted to real world situations than the RSB calculations. Such RSB calculations of 
fluctuations around the mean-field solution have been attempted for the random-manifold problem only in 
the case of 1-step non-marginal RSB (with disappointing result [58]) and offer, in full generality, extreme 
complications, as is illustrated by several studies for spin glasses [59-66]. 

The next challenge is thus to extend the FRG in a large-TV expansion beyond the dominant order (TV = 
oo). This is the aim of the present paper. Since this is a complicated calculation, and involves developing 
new methods which are of interest by themselves, this paper is restricted to the calculation of the effective 
action and derivation of the /^-function to order 1/iV. This is performed at T = and at finite temperature. 
The analysis of the resulting FRG flow, comparison with other studies, and applications, e.g. to the strong- 
coupling phase of the KPZ equation is involved and is the subject of a forthcoming publication. 

The outline of the paper is as follows. In section 2 we give the general formulation of the 1/JV expansion 
for the effective action of the random manifold. Details and generalizations are given in appendices A, B, 
and C. In section 3 we summarize the main results for TV = oo. Section 4 explains the derivation of the 
1/iV correction by a graphical method, which introduces a new type of diagrammatics. Section 5 explains 
the principle of a second and complementary method based on the algebra of 4-replica tensors. Section 
6 contains the full result for the effective action to order 1/TV, first expressed in bare parameters, then as 
a function of the renormalized dimensionless disorder. This allows, in section 7, for a derivation of the 
/3-function at T = 0. The structure of the finite-T /3-function is indicated, and details given in appendix H. 
A fool-proof diagrammatic version for finite temperature is given in appendix D. More details on the two 
main methods are given respectively in appendix G (for the diagrammatic method, including an alternative 
derivation of the T = /3-function) and in appendix F (for the algebraic method). Appendix I contains a 
list of all integrals. A table summarizing the notation is found in appendix J. 

2 1/TV expansion of the effective action: General formula 

We start from the partition function of an interface Z v = f P[ii] e~ Hv ^ u ^ T in a given sample, with energy 

Wv[«] = 7>(q 2 + m 2 )u(-q) ■ u(q) + V(x,u(x)) , 



(2.1) 



where J = J -^S^, J x = J d d a; and u ■ v — ^2 i=1 u l v\ The O(N) indices will be specified only when 
strictly necessary, and below additional replica indices for the replicated field u l a will be introduced, a = 



1, . . . , n. The small confining mass m provides a scale. To obtain a non-trivial large- N limit, one defines the 
rescaled field v = u/\fN and chooses the distribution of the random potential to be rotationally invariant, 
e.g. its second cumulant as 



V(x, u)V(x>, u>) = R(u - u')8 xx , = NB((v - v'Y)8 xx , (2.2) 

in terms of a function B(z). Higher connected cumulants are scaled as 

V(xi,ui) ...V(x p ,Up) = NS^^.^S^tvi^.^Vp) , (2.3) 

with 5 Xl: ... >Xp : = nL2 $ d ( x i - x i)- 

Physical observables can be obtained for any N from the replicated action at n = with a source 
J = y/Nj as 

Z[J]= f V[u]V\x}V[X}e- N ^ u ' x ' x ' j] (2.4) 

S[u,x,Kj] = ^ (q 2 + m 2 )v a (-q) -v a (q) 

+ / U(x(x)) - -iXab(x) [Xab(x) - v a (x) ■ v b (x)} - j a (x) ■ v a (x) , (2.5) 

where the replica matrix field x( x ) — Xab( x ) has been introduced through a Lagrange multiplier matrix 
field Xab(x). Here and below summations over repeated replica (and O(N)) indices a, b = 1, . . . , n is 
implicit. The bare interaction matrix potential 



U (X) =~7^Y1 B (Xab) ~ 7^ Yl S (* ab > * bc > ^ + • • • (2 - 6) 

ab abc 

depends only on the matrix 

Xab ■ Xaa + Xbb Xab Xba \^- ' ) 

and has a cumulant expansion in terms of sums with higher numbers of replicas. 

The effective action functional T[u] is defined as the Legendre transform of W{J] = \nZ[J] and 
satisfies 

I>] + W[J] = f J(x) ■ u(x) . (2.8) 

Since T[u] defines the renormalized vertices, its zero-momentum limit defines the renormalized disorder, 
the quantity on which we focus here. Thus we only need the result (per unit volume) for a uniform config- 
uration of the replica field u a (x) = u a = \^Nv a , which takes the form: 

m := ^r(u) = ^m\l + U(vv) , (2.9) 

where vv stands for the matrix v a ■ v t- This defines the renormalized disorder potential U{vv) and, whenever 
it can be expanded, up to a constant, 

U(vv) = -± J2 Btfb) -^lE § (&, "I, v 2 J + -... (2.10) 

ab abc 



It defines the renormalized cumulants B(z), S(. . .) etc.. Here and in the following we denote 

Vab-=V a -V b . (2.11) 

We aim at calculating the effective action up to terms of order 0(l/N 2 ), i.e. the first two terms in the 
expansion: 

U(vv) = U°(vv) + ^U\vv) + 0{l/N 2 ) . (2.12) 

Details of the calculation, as well as expressions for non-uniform fields are given in Appendix C. For the 
leading term we find, from a saddle-point evaluation [8]: 

U°(vv) = U(xv) + 7i E "TT^ 111 l-2Td x U( Xv )T +1 (2.13) 

n=l 
In '- = A(Jfc 2 +m2)n' (2 - U) 

The trace acts on replica matrices and \v satisfies the self-consistent equation 

„ oo 

Xt = v a v b + T Gt(k)=v a v b + TI 1 5 ab + Tj2 I n+i(i-2Td x U( Xv )} n )a b (2.15) 

^ k 71=1 

G v (k) = [(k 2 + m 2 )5 + 2Td x U( X v)] "' • (2.16) 



Note that for d < 2 no UV cutoff is necessary (apart for a constant term in the free energy), while for 
2 < d < 4 an UV cutoff is necessary (and implicit in the following) only for I\ x . 

One also finds a compact and very useful self-consistent equation for the derivative of the zero-th order 
potential: 

d ab U°(w)=d ab U( Xv ) 

= d ab u(vv + TI l 5 + TY j I n+1 (-2TdU\vv)) n \ . (2.17) 

Everywhere we denote by d ab U((f)) := d^ a JJ{4>) the simple derivative of the function U((j>) with respect 
to its matrix argument 4> ab . (Note that d ab U{vv) is a first derivative of U(vv) with respect to the matrix 
element v a ■ v b .) 

Next, from calculations of the fluctuations around the saddle point, one obtains the \/N correction, 
which can be expressed in terms of the zero-th order quantities as: 

U\vv) = \ |tr (In [5 ac 5 M + 2Td Xab d Xef U( Xv )(TU e J' cd (q) + v e Gl c (q)v d + v f G e v d (q)v c )]) (2.18) 

K f ' cd (<l)= I Gf{k)Gi c {q - k) (2.19) 

Jk 

G v (k) = [{k 2 + m 2 )5 + 2Td vv U°(vv)}- 1 , (2.20) 



To obtain a correct continuum limit, T therefore should be scaled as T = T/A , when A is taken to infinity 
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Figure 1: Top: typical T = contribution to B(v a b)- Bottom: self-consistent equation at leading order for B'(v^ b ) — B'(x a b)- 
The wiggly line denotes a derivative, and is combinatorially equivalent to choosing one B. At finite T one can attach an 
additional arbitrary number of tadpoles to any B. Also note that no loop made out of 3 propagators appears: this would be a 
contribution to the third cumulant (3-replica term), not calculated here; it is given in [8]. 

where here the trace acts in the space of replica pairs, i.e. tr(M) = J2 a b M a b,ab- Note that U°(vv) can also 
be replaced by the full U(vv) in the expression of U x {vv) with the same accuracy (i.e. at leading order in 
l/N). 

The saddle-point equation (2.17) for the zero-th order and the result for the 1/iV-correction (2.18) are 
still formal as they encode the full renormalized disorder distribution. To yield the renormalized disorder 
cumulants via (2. 10), (2. 12) they must be expanded in the number of replica sums, i.e. in cumulants. In the 
following Section 3, we recall the results for N = oo, and proceed with the non-trivial evaluation of (2.18) 
via a graphical method in Section 4, and via an algebraic method in Section 5. 



3 Review of the results for N 



oo 



In this section we review the main results at N = oo. Details can be found in [6]. 

3.1 Self-consistent equation at N = oo 

We start by recalling the cumulant expansion and only derive the result for the second cumulant. Higher 
cumulants are given in [6]. One studies a bare model, where only the second cumulant is non-zero: 



u(x) 



2T 2 



J^Z2 B (Xab) 



(3.1) 



ah 



and calculates the renormalized disorder (2.10). (We will drop the index zero on the cumulant functions, 
indicating the leading order). The self-consistent equation (2.17) can be expanded in sums with increasing 
numbers of replicas. We only need: 



-2TdU°(vv) 
-2TdU°(vv)) 2 



7f;(Sa b ^ & ac-B' ab ) + ... 



T 



my [&<*?] KeKf ~ B 'ab z2( B 'af + B 'bf) + z2 B ac B cb 



J-2 



ef 



(3.2) 
. . . , (3.3) 



where B' ab = B'(v^ b ) (recall v ah := (v a — Vb) ). Here and below the dropped terms contain sums with 
too many replicas to contribute to the final result for the self-consistent equation of the second cumulant 



(2-replica term). One thus has: 



+ 2TJ 1 (1 - S ab ) + 4/ 2 [- J2(K C + B' bc ) ~ Ma + B'J + B' ab - S ab £ B' a 



X v =v a 



The self-consistent equation becomes: 



+ 



(3.4) 



B'(v 2 ab ) - S ab J2 B'(vl) + 3-replica terms = B'{xf) - 5 ab £ B'(xT) (3.5) 



and can be solved by appropriate Taylor expansion of the r.h.s.. It is solved for a ^ b: 



B' (v 2 ab ) + i J2 S'<*9 + E ••• + ••• = B' (xt) ■ 



T 

9 gh 

It is then easy to see that the second cumulant satisfies a closed equation at any T, 



(3.6) 



B'(v 2 ab ) = B' v 2 ab + TTh + 4h(B>(v 2 ab ) - B'(0)) , 



(3.7) 



with no other contributions from higher cumulants at any T. A more detailed derivation is given in [6]. 

3.2 Derivation of the FRG equation at N = oo 

From the previous section the renormalized second cumulant of the disorder B'(x) satisfies the self- 
consistent equation 

B'(x) = B' (x + 2Th + ±h{B'{x) - B'(0)j) . (3.8) 



It implies 
as well as 



B'(0) = B'(2TI 1 ] 



B"{x) = B"{x + 1TI X + 4/ 2 (S'(x) - B'(0)))[1 + 4I 2 B"(x)] . 
We now derive the corresponding exact FRG equation. Taking the derivative md m gives: 

md m B'(x) = B"(x + 2TI X + AI 2 (B'(x) - B'(0))) 

2Tmd m h + 4(md m I 2 )(B'(x) - B'(0)) + M 2 md m B\x) - 4I 2 md m B'(0) 
B"(x) 



(3.9) 
(3.10) 



x 



l + 4I 2 B"(x) 

x[2md m TI 1 + A(md m I 2 )(B'(x) - &(0)) + AI 2 md m B\x) - U 2 md m & '(0)] . (3.11) 



This yields: 

md m B'{x) = B"(x)[2md m Th + A(md m I 2 ){B \x) - B'(0)) - AI 2 md m B' 

Thus one has also: 

md m B'(0) = ^-22 2md m (TI 1 ) . 
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(3.12) 



(3.13) 



■O <!> <0 {4) o 



Figure 2: A typical "tree-like" diagram contributing to the leading order in the 1 /TV-expansion of the renormalized potential V 
(effective action). The numbers 1, 2 and 3 depict special features discussed in the main text. There is a factor 1/7V per dotted 
line and a factor of N per "small loop" of propagators (solid lines), thus the overall factor is N as it should from (4.1). This 
graph does not contain any "big loop" (see text). 



Hence one gets finally 



m^-B'ix) = B"(x) 
dm 



2( m AtM L + M m JLl 2 )(B'( x ) - B'(0)) 

K dm J 1 + 4I 2 B"(0) dm 2JK K J K " 



(3.14) 



which can also be integrated once over x. As emphasized in Ref. [6] it is exact at N = 00 for any d and 
correctly matches the 1-loop FRG equation obtained by Balents and Fisher for any N but only to 0(e), 
e = 4 — d. It can be solved directly, or equivalently the self-consistent equation (3.8) can be inverted. 
The corresponding solutions for various models are discussed in [6] and compared with the Mezard-Parisi 
solution [11] obtained in a rather different manner through a replica- symmetry-breaking saddle point. 

Before discussing specific models, we now turn to the evaluation of the effective action of the FRG and 
the /3-function to the next order in 1/iV. 

4 Corrections in 1/JY, via the graphical method 

In this section, we present a graphical method to calculate the corrections to B at order l/N. An alge- 
braic method is presented in the next section 5. Both methods are completely independent, since they use 
orthogonal ideas. They were performed independently by the authors, each on a different continent. The 
agreement on the final result gives some confidence that it is free from calculational errors. 

4.1 General considerations for a scalar field theory 

Let us start with some general considerations about which graphs contribute at a given order in l/N. For 
that purpose, we consider a general (pure, no disorder) scalar field theory, with a iV-component field u l (x), 
i — 1...N and interaction 



S** = NjV\z r ) , (4.1) 

where the wM-correlations of the free theory (V = 0) are given by 

(u i (x)vP(y)) = S i 'C(x-y), (4.2) 



with i,j = 1 . . . N and C(x — y) independent of N is denoted by a solid line in figure 2. Graphically, we 
can denote this by 

^ = h I + I ) { + W 2 



+ ... , (4.3) 



where all coefficients bi are of order 1; there is thus a factor of 1/JV per dotted line. 

The renormalized potential V (effective action) is given by the sum of all 1 -particle-irreducible di- 
agrams. They must thus contain "small loops". To leading order in 1/JV the following diagrams are 
possible: 

• Tadpoles, contracting any u 2 with itself only. The factor of N from the "small loop" (the ^2 i=1 over 
the number of components) compensates the factor of TV from the argument of V(u 2 /N). (See 1 on 
figure 2.) 

• Closed "small loops" with 2, 3, 4 or more vertices, as denoted by the same number on figure 2: 
Adding one more vertex gives a factor of N, see Eq. (4. 1); this is compensated by using an additional 

u 2 /N. 

• Note that all vertices u 2 , (u 2 ) 2 , a.s.o. contribute equivalently. (See 5, 6 on figure 2.) 

Note that all these diagrams are "tree-like" diagrams, where the branches (made of "small loops") are 
made out of the diagrams through which no total momentum is running. We will call them "tree-like" in 
the following, to distinguish them from normal trees. They are resummed by the saddle-point equations 
(see appendix B) or from graphical inspection, as done here: First of all, insertions of V into a line of 
propagators act like a mass, leading to the replacement of l/(k 2 + m 2 ) by 

G(k) := l . (4.4) 

k 2 + m 2 + 2V'(j r ) 

The effective potential V is obtained from 

V' (f ) = V ( X ) (4-5) 

JV J k JV J k k 2 + m 2 + 2V I ^ 



N 

Note that the derivatives are graphically understood as follows: Choosing one vertex (derivative!) in the 
effective potential V is equivalent to having a bare vertex with the same derivative taken (thus V) and 
attaching to it loops made out of correlation-functions. Attaching any number of such loops to V, amounts 
to shifting its argument, as can be seen from Taylor-expansion. In these loops, again derivatives (one needs 
u 2 to attach the loop) of the effective potential are inserted. The latter can thus be written as V with shifted 
argument to account for more things to be attached to this V or equivalently using (4.5) to a V'. This result 
coincides with (B.22). 

Diagrams at next order 1/JV contain exactly one "big loop", see figure 3. Take a "tree" (as e.g. the 
object on Fig. 2) and glue it together to form a "big loop" by identifying two vertices; this does not change 
the factor of N from the loops, but one looses one factor of JV from the missing vertex. Note that also the 
"minimal loop" marked on figure 3 belongs to the same class, even though it looks different. The "big loop" 
demands to sum a series of diagrams at non- vanishing momentum, and then to carry through the integration 
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minimal 
loop 



Figure 3: Typical graphs which contribute to the renormalized potential V to subleading order in l/N. The "big loop" (a loop 
made of loops, see main text) accounts for a factor of l/N; the same is true for the "minimal loop" on the right, which is a "big 
loop" in disguise (as explained in the main text.) The given diagram is thus of order l/N 2 

over momenta. However the simplification remains that any added "tree-like" branches are resummed by 
replacing the argument of V from u 2 /N to x- 

Another feature arises at subdominant order: In the "big" loop, one may pick any given small loop (i.e. 
the loop made out of correlation-functions) and replace one of the two correlation-functions by u 2 /N: This 
means that the corresponding fields u did not get contracted. Since the remaining correlation-functions 
force their indices to be equal, this gives a factor of u 2 = N x [u 2 /N), thus contributes the same factor 
of N. Note that these diagrams do not contribute to the effective action at leading order, which is treelike, 
since the resulting diagrams would be 1 -particle reducible. 

To resum the order 1/iV-diagrams one has to sum over loops of all sizes. The result is 



SV[% 



1 Al 



£ 



2N^n dl> 

71=1 r 



2V"(x)(h(p) + 2fG(p) 



= -1 |ln [l + 2V"( X ) (hip) + 2 f G(p) 
Hp)= ( G{k + P )G{k). 



(4.7) 
(4.8) 



This can be compared to the results of appendix B, and more specifically to formula (B.39). 



4.2 Elastic manifolds in disorder: General considerations, building blocks 

Let us start the treatment of the disordered model with some general considerations. First to organize the 
1/iV-expansion, one may still use the diagrammatics of the previous section, which shows the 0{N) -index 
content. The same diagrams still exist, but they now also have a complicated replica content. The replica 
content can be explicated by using "splitted vertices" instead of the unsplitted ones of the previous section. 
The corresponding replica diagrammatics, which shows the replica structure only was explained in details 
in [7]. This can be drawn as 

Y,B{{v a -v h ) 2 ) = a *- b , (4.9) 

ab 

where a dashed line connects the two dots, standing for replicas a and b. In order to avoid confusion, note 
that this dashed line is different from the dotted line used in figures 2 and 3 as well as equation (4.3) to 
show the O(N) -structure. 







Figure 4: Example of a big loop made out of small loops, in splitted replica notation (left). It contributes to the 2-replica part (i.e. 
the disorder) at order 1/iV. Note that there are two constraints more than needed to have a 2-replica term. Thus two redundant 
constraints can be "wasted", by cutting each of the solid lines exactly once, either by inserting a B' or leaving a uu uncontracted, 
as is done on the right. 

Below we introduce a third diagrammatics which allows to track both the replica and the O(N) indices, 
not an easy task. Before doing so let us explain a few points. 

Since we are only interested in the corrections to the 2-replica part B ab of the effective action, there are 
many O(N) diagrams which do not contribute. 

At dominant order (figure 1) small loops with three vertices (see 3 in figure 2), or more, do not contribute 
to B ab , but to the third cumulant, or higher, as can be seen from 



(4.10) 



which is a diagram showing only replica indices, N indices being implicit. Note that solid lines iden- 
tify replica indices. This is why in figure 1 only chain-diagrams and tadpoles (the latter are omitted for 
simplicity of presentation) appear, rendering calculations appreciably simpler. 

At subdominant order in l/N, only diagrams with one "big" loop (general feature of the order l/N 
discussed in the last section) made of any number of "small" 2-loops and only exactly two "small" 3- 
loops can contribute to the renormalized second cumulant. Each small 3-loop can also be replaced by an 
uncontracted uu. To understand this, consider the simplest "big" loop, i.e. the railroad diagram, which is 
drawn on the left of figure 4 in splitted replica notation. It contains exactly two closed propagator lines, 
which over-constrain the replicas to be equal. (These are the inner and outer solid lines on the left of figure 
4). These over-constraints can be relaxed, by cutting each line exactly once, in order not to get a higher 
replica term, as is illustrated on the right of figure 4. This "cutting" is possible by either inserting into a 
propagator a vertex (which contains two replicas, thus is not "replica-conserving") or by leaving one uu 
uncontracted. This is the basic principle, whose careful exploration leads to all of the diagrams at order 
l/N, as we will discuss now. 

In order to do so, we have to introduce a new powerful graphical notation: 

■ =B(x) (4.11) 

- a b =B\x ab ){u a -u b f (4.12) 

\ ^ a h =B"{ Xab ) l -[{u a -u b f]\ (4.13) 

where lines departing in the same direction belong to the same vector-index, and a continuing line repre- 
sents the same replica. 

We also use the following short-hand notation 

B' ab :=B'(xa b ) (4.14) 
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and similar formulas for the higher derivatives B" ah := B"(xa b ), etc. Another frequently used shorthand is 

B' :=B' aa , B»:=Bl, etc.. (4.15) 

In order to be able to resum the "big" loop, we will now introduce some building blocks. Since there can 
be any number of 2-loops in the "big" loop one needs to define the resummed chain 

;o; :=: = ;+:o;+:oo:+:<xx>:+-- <4.i« 

We have with momentum p running through the diagram 

^ a b O a b =4h(p)B"(xa b ) a b ^ a b (4.17) 

with 



h ^ J k (k + p/2) 2 + m 2 (k - p/2) 2 + m 2 ' (4 ' 18) 

Let us introduce a compact notation for the integrals, summarized in appendix I. (All our notations and 
important formulas are also summarized in a table in appendix J.) 

f 1 1 

Jaf){q) = J q a( ] ■= J^ ( (A; + g / 2 )2 + m 2y (77. _ q / 2 y + m 2y ( 4 -!9) 



1 1 

k {(k + p/2) 2 + m 2 ) 2 (A: -p/2) 2 + m 2 



hip) := J,M = I ,„ , —, , ^7r—j7^—- 2 (4-20) 



hip) ■■= h lip) = f o o • (4.21) 

V ; ' K ' Jk iik + p/2) 2 + m 2 ) 2 Hk - p/2) 2 + m 2 ) 2 



Thus 



<x> : = MP) 5"(^)] 2 ; - : (4.22) 



-^ ; <xx> ; = [4/ 2 (p) 5"(^)] 3 ; - ; , (4.23) 

and so on. These chain-like diagrams form a geometric series, which is resummed as 

:o: = Lu^B-i^) a ^>'-\ [{Ua - Ub)2] Hab{p) [{Ua - Ub)2] ■ (4 - 24) 

We have introduced H ab (p), the "effective" B" ah after resummation 

H ab ( P ) := : A B '}fi ( - s , (4-25) 

1 - AI 2 ip) B"ixab) 

which we equivalently can express at leading order in I/TV through B"(x ab ) = B" b /[1 + 4I 2 B" b ] as 

H ab (p) := ^ ^- . (4.26) 

l + 4[/ 2 -/2(p)]^ 
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B(Xab) 
B'(Xab) (Ua - Ub) 



B"( X abh[(u a -U b ) 2 } 



21 2 



SabCip) 



>ab 



p2 _j_ m 2 



''ab 



a , — N a 
b A/ h 



H ab (p) 



B"( X 



ab 



4I 2 (p)B"(xo 



a ^^ b 



8 ab I 2 (p)(l - ±I 2 B>> a 



a ^-^ a 



I 2 (p)(l + 2I 2 B': a )(l-AI 2 B'; a 



Figure 5: Building blocks of the perturbation theory. See main text. The last two blocks only appear at finite temperature. 



We also define H v (p) as 



H v (p) 



B"(x(v)) 



B"(v 2 



1 - AI 2 (p)B"(x(v)) 1 + 4 [I 2 - hip)} B"(v 2 ) 



(4.27) 



valid again at leading order in l/N. Note that the denominator in (4.27) reflects the renormalization of B: 
The divergent integral I 2 (p) does not appear alone, but together with its counter-term, the integral I 2 (p) 
subtracted at p = 0. It turns out that at zero temperature the above are the only building blocks needed. 
However at T > one needs two more building blocks which are quite non-trivial. As was shown in [7] 
non-zero temperature diagrams contain at least one replica "sloop"; these are exactly the over-constraining 
lines discussed above and on figure 4. There is one factor of T for each such "sloop". The additional 
building blocks thus contain sloops: One sloop at order T, and 2 sloops (an example is on the left of figure 
4) at order T 2 . Higher orders in T are only possible at order l/N 2 , or higher. 

To explain the construction of these additional building blocks, which is subtle, one goes back to the 
diagrammatics showing only replica indices. 

The first building block is the "moon-diagram". This is the sum over all diagrams, which at both ends 
have lines joining only one of the both replicas, and which enforce the joined replicas to be equal: 



^=^ b 




(4.28) 



Note that we construct the chain from left to right. Otherwise the graphical representation of the perturba- 
tion expansion is not unambiguous, as can be seen from the following example 

(4.29) 
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Drawing first the two left-most lines, the two right-most ones can no longer be added, since in the middle, 
there is R"(u a — u a ) = R"(0), which does not depend on the field. Conversely, if we decided to first draw 
the two rightmost lines, then the two leftmost could be drawn. Consequently the diagrams to be drawn 
in (4.28) would be different (actually they would be nothing but the diagrams mirrored such that their 
left and right ends are exchanged), even though the final result would be the same. This phenomenon is 
detailed (for a different diagram) in appendix E. Note that there is nor a contradiction, nor an inconsistency 
of the approach. It merely means, that when using these kind of rules, which have the advantage of 
simplifying calculations importantly, one has to order the contractions. An approach, which does not have 
this deficiency, but is very complicated, is explained in appendix D. We will use it there to recalculate 
diagram (4.28). 
We claim that 



^=\ 



4/ 2 (pK a 



(4.30) 



The second identity is trivial perturbation theory. The non-trivial statement is the first identity. To prove it, 
we remark that starting with (recall we construct from left to right) ; ^^^ no further contraction can be 
made. Therefore, we have to start with ^ — ^>» • At chain-length two, the only possibility is 

(4.31) 



At chain-length three, there are two and only two possible prolongations, which have no additional free 
replica-indices: 



0. 



(4.32) 



These diagrams cancel. The same is true for longer chains, since at any intermediate position (i.e. not the 
first and not the last lines), there is always the combination 



(4.33) 



which when closed at the right end with ; ^^^ ; cancel. This completes the proof. 

The last diagram which we need is the "half-moon-diagram" /^\ , which is similar to the moon- 
diagram, but does not enforce the replicas at its ends to be equal. However it will always be evaluated at 
coinciding replicas. (It thus contains ^=\ as a subset.) We claim that 



a ^^ a 



{l + 2I 2 (p)Bl)(l-AI 2 (p)Bl). 



(4.34) 



This is proven by first remarking that all other diagrams can be generated from those: Add left of a £ 
or a !v~^x! the combination 



they cancel pairwise. The same is true for chains made out of the combination ■ ■ + 

note that one cannot insert more Jv^n; , since they would lead to higher replica terms. 
All rules and building blocks are collected on figure 5. 
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4.3 Zero temperature (T = 0) 



We start our discussion with zero temperature, T = 0. We have to construct all diagrams with the topology 
of a loop in the large- N limit. Note that e.g. the diagram (4.17) counts as a line in this construction. The 
building blocks are given on figure 5. At zero temperature, one needs all possible constraints on the sum of 
replicas. This means that one can not use ^=\ or /=x , which both contain one non-replica conserving 
line !<~>I a t zero temperature. At finite temperature, one can use one at order T and two at order T 2 . 

At T = 0, we find the following diagrams. Note that the notation is such that crossing lines do not 
intersect. All the diagrams correct the effective action B without any further combinatorial factor. 




2 12 



8 / C(p) 2 [H ab (p)v 2 ab ] 



(4.36) 



64 f I 3 (p)C(p)H ab (p) 2 v 2 ab (B' ab - B' ) (4.37) 

Jv 



128 fl 3 (p) 2 H ab (p) 2 (B' ab -B t 
Jp 

8 [c(p) 2 v 2 ab H ab (p)(B> ab -B> 



o) 



16/ H ab (p)(B' ab -B' yi 4 (p) 

J k Jp 



(4.38) 
(4.39) 

(4.40) 



The combinatorial factors can, and have been checked by straightforward calculating the diagrams with 
H ab replaced by B" b , both by hand and computer-algebraically. 

The idea of how to construct these diagrams is straightforward: We start by a closed chain of 2-loops, 
which disregarding the O(N) -structure has been drawn on the left of figure 4. Then one has to cut each 
line exactly once. These cuts can either be done at different positions in the "big" loop (diagrams (4.36) to 
(4.38)), or at the same position (diagrams (4.39) and (4.40)). Then there is the possibility to either insert into 
a propagator a vertex (these are the terms proportional to B' ab — B' ) or not to contract two fields (the terms 
proportional to v 2 b ). The left and right diagrams in each equation are distinguished by their twist: The left 
one is the untwisted one (as the one drawn on the left of figure 4), the right one the twisted one, obtained 
by cutting both lines between two neighboring vertices, and reglueing them together with the two lines 
exchanged (this gives one single propagator line running twice around.) Note that the twisted diagrams do 
not appear in the final result, since we suppose analyticity for B, such that e.g. \\m. a ^ b H ab v 2 b = 0. 2 

As discussed above other imaginable contributions are 3-replica terms, where one is not using the 
maximal number of possible constraints on the number of free replica-sums (one line cut twice instead of 
each line cut once): 




3-replica term 



(4.41) 



2 Note that this construction suggests how to construct additional "anomalous" terms, known e.g. to be necessary at 2-loop 
order [7,50]. We do not present them here, but relegate their discussion to a subsequent publication. 
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3 -replica term 



We note contributions which vanish identically for completeness: 




(4.42) 




0. 



(4.43) 

(4.44) 



4.4 Corrections at order T 

We remark that adding an additional line to any object already constructed at order 1/iV results into a 
diagram of higher topology in the large iV-limit. This means that the diagram does not contribute, a 
statement which remains true to any order in perturbation theory. The only remaining possibility is to 
proceed as before, but constraining replica-indices to be the same by more than one propagator (line). 
Examples are ^=\ and ^=\ . Another example would be a "circular railroad diagram", see left of figure 
4 and equation (4.50) below. Since we had two lines to "waste", it means that there will be a term of order 
T and T 2 . (As a side-remark, we note that at order l/N 2 , there will be terms of up to order T 4 , since for 
the leading term one has to cut up to 4 lines, a.s.o..) 
We now give the order T-contributions: 



AT 



v 2 ab Hab{p) 



jr 



m 




16T I 3 (p)H ab (p)[B' ab -B' ]. 



(4.45) 



(4.46) 



These are the contributions, where one replica-line has been cut by the insertion of either v 2 or B', whereas 
the other one (on top of the diagrams) is redundant. 

We can also use a double (redundant) line, using the moon diagram ^=x . Starting from Eq. (4.45), 
and inserting it into the bubble line, we find 



Q^) = 8T J C(p)vl b H ab (p) 2 I 2 (p) [1 - AI 2 (p)Bl 



(4.47) 



Starting the same procedure from (4.46) gives 



<C^C^) = 32T / H a b (p)\B' ah - B' aa )I 2 (p)h(p) [1 - AI 2 (p)Bl 
Note that there is also a "twisted" version of (4.46) and (4.48), which add up to 



(4.48) 





16T H aa (p)B' ab h(p) [l + 2I 2 (p)Bl 



(4.49) 
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4.5 Corrections at order T 2 

We continue with diagrams at order T 2 . There is one diagram, which does not necessitate any cut in a 
replica-conserving line, thus has two redundancies. It is the "railroad-diagram", i.e. a closed chain, with 
at least one vertex. Note that contrary to what one might expect, for one vertex this is not a product of 
two tadpoles summed at leading order, even though it looks alike. However, the structure in vector-indices 
makes it a loop in N, i.e. a subdominant term. Since no vertex is marked, the sum is not a geometric series, 
but a logarithm: 

= g ,£ im-wir (450) 

This can be written either as a function of the bare disorder, or using B"(xab) = B" b /(1 + AI 2 B'^ h ) (valid 
at all T) as a function of the renormalized disorder: 






l-4J 2 (p)£"(x t 

l + AB"(v)[I 2 -I 2 (j>)] 
l+4B"(v)Io 



(4.51) 




Note that this is the term where a double, completely replica-conserving line goes around. 

More diagrams are possible, with one and two defects, i.e. replica-line cuttings, equivalent to inserting 
/^=x or /=\ . Using only one cut, both outer indices are forced to be equal, and we have to insert /=\ , 
calculated in (4.34): 

IT 2 fl 2 (p)H ab (p) [1 + 2I 2 (p)Bl) [1 - AI 2 (p)Bl] . (4.52) 

The overall prefactor of Eqs. (4.51) and (4.52) is such that the term linear in B" ah comes with a factor of 2, 
and they both add up to a factor of 4, which can be checked with a simple 1-loop calculation. 

The remaining term is obtained by cutting two-replica lines, using the replica-conserving moon ^=\ . 
This gives the contribution 

= IT 2 I [I 2 (p)H ab (p)(l - AI 2 (p)Bl)] 2 . (4.53) 

Jp 

There is an additional anomalous term. It is nothing but a 1-loop diagram, of the form ; ^^^ ; , where 
the right-most vertex is the sum of all diagrams at order T 2 , as given by (4.50) to (4.53), evaluated at 
coinciding replicas: 

A i° lftT 2, R , f h(p)[l + 2I 2 (p)Bl]BZ 
A =- 16ThBab J P l-4/ 2 (p)^ 

= -lQT 2 I 2 B' ab fl 2 (p)[l + QI 2 (p)H aa (p)]BZ . (4.54) 

Jp 

This can be reexpressed as a function of B, see section 6.2. One might suspect that longer chains can be 
constructed to connect the B' ab with the derivative of the effective action at order T 2 , taken at coinciding 
indices. With the same arguments as already made a couple of times above, the insertions of ■ ■ + 
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; ^^^ ; pairwise cancel. Also note that one could of course draw all these diagrams; there is however no 
clear advantage of doing so. 

The alert reader will also wonder why we have not mentioned any such term at order T° or T. It turns 
out that the effective action at order T°, when derived once and taken at coinciding arguments, actually 
vanishes (supposing analyticity!). At order T, the term in question is nothing but (4.49). 

One caveat is in order: Even though this procedure is simple and elegant, there are many hidden traps. 
It is therefore good to check this calculation by an explicit loop expansion, using the excluded replica 
formalism. This has been done up to 8-loop order, and relies heavily on computer algebraic support. The 
procedure can also be formalized, leading to an additional more rigorous but somehow elaborate approach, 
the "Excluded Replica Formalism", which is presented in appendix D. 

5 Corrections at order l/JV, via the algebraic method 

The calculation of the renormalized disorder to next order requires the calculation of the following trace in 
the space of four replica matrices, from (2.18) and (2.19): 

U^vv) = - tr[\n(S ac S bd + M ab:Cd )] (5.1) 

M ab ,cd = M ab , ef N eUd (q) (5.2) 

_ M ab , cd = (2Td Xab d Xcd U( X ))\ x=Xv (5.3) 

N eUd {q) = TUf cd (q) + G f v c (q)v e ■ v d + G e v d (q)v f ■ v c . (5.4) 

n^' cd (g) has been defined in (2.19). Expression (5.1) can be computed by systematically expanding in 
sums over increasing numbers of replicas, as was done at dominant order. However the calculation is 
considerably more tedious in the present case. We give the main features here and relegate details to 
Appendix F. 

To obtain the correction to the second cumulant of the disorder, we will only need the two-replica part 
of this function, which we denote P2U 1 (vv) (P n X denotes the part of an expression X which contains 
exactly n free sums over replicas). Since the trace already involves at least one replica sum, we can and 
will truncate all expressions given below by discarding all terms with two or more replica sums. This can 
be checked systematically and originates from the fact that once a replica sum appears in an expansion, 
it can never disappear later on. As a result, we find that third and higher cumulants do not appear in the 
correction to the second cumulant, as they involve higher replica sums 3 . 

The expansion of the matrices M, N and M. are computed in Appendix F. It is crucial to write explicitly 
all Kronecker-delta functions. The matrix M. is found to have the form: 

M abtCd = m ab -(5 ad 5 bc + 5 bd 5 ac ) + (1 + m ab )M a b,cd , (5.5) 

where m ab is symmetric in a, b and A4 is symmetric in a, b and symmetric in c, d (but not necessarily in 
exchange of (a, b) with (c, d)) and can thus be parameterized as 

M-abfid = SabcdXa + ^abcUad + ^abdUac + $acdZ a b + $bcd z ba + 8 ac tabd + $adtabc + 8 bc t bad + 5 bd t bac 

+$ab$cdU a c + 5 ab V acd + 5 cd W abc + g a bcd , (5.6) 



3 This formal expansion in replica sums assumes some analyticity property in a way which should be analyzed later. 

17 



where all Kronecker-delta's have been written explicitly. Also note that v acd is symmetric in c and d; w a bc 
is symmetric in a and b and g abcd is symmetric in a, b and symmetric in c, d, whereas all others have no such 
symmetry. These matrices form a closed algebra which is studied in Appendix F. Unfortunately this algebra 
is rather large, even though it is the smallest algebra sufficient for the present calculation. Note that we 
have explicitly separated the part proportional to the (symmetrized) identity l^™ d = l/2(5 ac 5 bd + 5 ad 8 bc ). 
The first preliminary step in the trace log calculation is to prove 4 that 

i i 

tr [\n(5 ac 5 bd + M ab ,cd)] = ^ Yl ln ( l + mab "> + 2 S ln ^ + maa "> + tr [ ln ( 5ac5bd + M ^,cd})] , (5.7) 

ab a 

where in the first two terms the In simply acts on numbers. This formula is valid 5 for a matrix M. symmetric 
in a, b and symmetric in c, d. The last trace log is equal to its usual series expression XL>i ^ — {M Yob ab- 
It does not seem possible to express the trace log in general but here it is possible to expand it sys- 
tematically in the number of replica sums. Let us sketch the method that we found most convenient. We 
write 

r 1 c ]\ 

tr[ln(l + M)} = - / — trM A (5.8) 

Jo * 

M x = (1 + XM)- 1 - 1 , (5.9) 

and we want to compute all terms of M A , namely: 

M ab,cd = $abcdX a + ^abcVad + ^abdVac + ^acdZ ab + hcdha + ^act abd + Sadtabc + ^bad + $bdhac 

+5ab5 c dU X ac + S ab v x cd + 5 cd wl bc + gl bcd . (5.10) 

We do this by using the algebra detailed in Appendix F to solve the following equation, equivalent to (5.9): 

XM + M A + A M M x = . (5.11) 

One projects onto each component x, y, z, . . . and onto terms with an increasing number of replica sums. 
What we want are x x , y x , z x , ... as a function of the known x, y,z, . . ., which parameterize M. (their 
expressions are given below). 
More specifically one writes: 

x a = P x a + P 1 x a + ... (5.12) 

x x a = P x x a + P 1 x x + ... (5.13) 

and similarly for all other components y, z, t, ... of the matrices M. and M x . Using the algebraic rules 
for the product of two matrices, it turns out that it is possible to solve (5.11) for all components of M x 
in an iterative manner. This is simplified since P^x a = (see below). First we determine all zero-sum 
components Pqx x , Poy x , . . . from the corresponding P x, P y, It can be done in the following order: 



4 Proof: Write M = A + (1 + A')M, where l abtCli = 5 ac S bd - Then tr[ln(l + M)\ = tr[ln(l + A + (1 + 
M)M)] = In (dct[l + A + (1 + M)M\) = In (dct[l + A])' + ln_(det[l + (1 + h)- l {t + A')M]). Further In (dct[l + A]) _= 
tr[ln(l + A)] = |E«i, ln ( 1 + m »t) + iE« ln ( 1 + m »»)- Since. Mab.cd = M ba .cd, the following identity holds: AM = A' M, 
such that In (det[l + (1 + A)" 1 (1 + A')M\) = In (dct[l + M\) = tr[ln(l + M)] = - E P >o(- 1 ) Ptr [ MP ] ■ q- e - d - 

5 Note that the naive identity tr[ln(<5 ac <5 bd + M a b,cd)] = ir[M^(S ac S bd + S ad S bc ) + M a b,cd)] is incorrect. 
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First compute PoVab together with Poy^, then similarly Poz x , Pox x , Pot x , Pov x , Potu x , Pog x - Second one 
can project (5.11) onto one-sum terms, and determine Pix x , Piy x etc. . . 
At then end we need: 



P 2 tr ln[l + At] 



A 



P 2 (trilP 



P 2 (tr M X ) = J2( P ^a + ZPWaa + 2Pl*&, + 2P^ aa + P lt £, + P lV ^ + P,W, 



J X ) 

J aaa) 



(5.14) 



ab 



,2Po^a66 + PoQabab) 



(5.15) 



The remaining integrations over A of each term are found to be of the form f dAA p /(l + 2y aa \) q = 
2 Pi(l + p, q, 2 + p, — 2y 00 )/(l + p) which for integer p and q can be expressed as rational fractions. 

The detailed calculation is however very tedious and has been performed using Mathematica. The 
result is given in the following and is found to agree exactly with the graphical method of the previous 
section. Here we give the zero-sum and one-sum components of the matrix M. needed for the calculation 
of P 2 tr ln(l + M). They are calculated from (5.1), see Appendix F, and read: 



P x a = , P u 
P y ad = h(q)Al a B 



-PoVa 



■ii 

ad 



P z ab = I 2 (q)A q ab B 

Po^abd — 
PoVacd = 



A-aa 



A q 

■ri-„ 



PoWacd 
2 







l-AB"(2TI 1 )I 2 (q) 
2 



--A q B" 



^ab^ab 



-A" l~R" 

rp aa I ad 

+ B 



' l-AB ab I 2 (q) 
(v b -v a )v d C(q)-2I 3 (q)(B' bd -B'J 

v d -v a )v c C(q)-2h(q)(B' dc -B' a 

v a )v d C(q)-2h(q)(B' cd -B' ad 



PoQabcd = T^C(q) 2 Al b B ab (v a - v b ) (& bc - B' ac )v d + (B' bd - B' ad )v c 



p2 



where we have defined B' 



ab 



^2' 



B\vl\) and 



h{q)A q ab B ab B' ad B' bc + B' bd B' ac - B' ad B' ac - B' bd B' b 



be 



B ai 

A v 



B"(Xa 



B 



ah 



B"\x 



ab) 

•jl b + 2TI l+ AI 2 {B' ab - l -{B' ac 



B'bb)) 



(5.16) 
(5.17) 

(5.18) 
(5.19) 

(5.20) 
(5.21) 

(5.22) 
(5.23) 



■ab 
v 



ab 
~(0) t 
Tab to 



with B aa = B"(2TI 1 ). All integrals I n and I n (q) are defined in section 4.2 and Appendix I and C(q) 
l/(q 2 + m 2 ). We recall that B is the bare second cumulant, B the renormalized one and satisfies B' a 
B'(j(v b ) at dominant order, which is sufficient for our purpose. For convenience we use x 
denote the zero sum part in the decomposition Xv b = Xab + $abXa m the notations of [6]. 
For the 1-sum terms we need only the diagonal values: 

P 1 x a = -I 2 (q)A aa J2Xf 
f 



(5.24) 
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PlVa 



, \B"(2Th)Uq) + TB / "(2TI 1 )I 2 I 2 (q)] V B' 

?"(2T/ 1 )/ 2 (g)f L V 1J SW V 1J iW J ^ a/ 



T [1 - 45 



/< 



(5.25) 



/ 



P\Zaa = -P\U, 

R 

PlVaaa = 



T 



[1 -4B"(2T h )h( q )? [B "^ 2TlMq) + ^(ZTJOW?)] E K f (5.26) 

/4(g)E^/(^/-^'(0)) 2 

(5.27) 



T 2 l-4 J B"(2T/ 1 )/ 2 (g)L 



+C( q rY / B 'af(B' af -B'(0))(v a -v f )v a . 

/ 

In addition since t aaa = one has P\t aaa = 0; PiW aaa = since w a6c = and it turns out that P\g aaa a is 
not needed. 

Starting from these values one performs the algebra and obtains P 2 tr ln(l + M). To this one must add 
the "simple" part of P 2 tr In in (5.7) above which one expands as follows: 



Po 



\llj ix ln{l + m ^ + ^ E / tr Xn{l + 



m n 



ab J 1 



(5.28) 



1 
2 



E Hi-aXMq)) 



ab J( i 



Tj q l-4B"(2Th)I 2 (q) 



I 3 (q)B"(2Th) + TI 2 B'"(2TI 1 )I 2 (q) 



B[, 



af 1 



where m ab is computed in appendix F. 

To translate the results into terms of the renormalized disorder, we can perform the replacements 

X b = K b /(l + AI 2 B': b ) (5.29) 

B"(2Th) = £"(0)/(l + 4/ 2 ir (0)) (5.30) 

Bl = BZ/{l + AI 2 K b f (5.31) 

B"\2Th) = B'"(0)/(1 + AI 2 B"(0)) 3 , (5.32) 

since to the same accuracy we can use in all above expressions the dominant order or the exact one. 
The "simple" part of P 2 tr In thus gives: 



5 {simplc) B(x) 



™ fh(p)Ho(p)&(x) + T -{ - l - /in fl + ^(x)[/ a -/ a (p)] 
N J p W)ow K) N \ 2 J p y l + 4B»(x)I 2 

+8/ 2 [ I 2 (p)B'(x) = ^^ = 

J P ; (l + 4/ 2 fi"(0)) 2 (l+4(/ 2 -/ 2 (g)) J B"(0)) 



(5.33) 



where integrals are defined in Appendix I. 



6 Results for the effective action 

The 1/N correction to the 2-replica part of the effective action in terms of the bare disorder can be written 



as: 



P*uHw) 



^E^o 



2T 2 



ab) 



(6.1) 
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with v 2 b = (v a — Vb) 2 . It can be read off from sections 4.3, 4.4 and 4.5. We will give the expression 
for the correction to the renormalized second cumulant SB(x) both in terms of the bare disorder and the 
renormalized one. 



6.1 The effective action as a function of the bare disorder 



We find: 



SB(vv) 



N 



21 2 



8 / C(p) 2 [H v (p)v 2 ] 



+ 



T 

N 



+64 / I 3 (p)C(p)H v (pYv 2 [B\x v ) - B'(xo)} 
Jp 

+ 128 / h(p) 2 H v {p) 2 [B\x v ) - B'(x )] 2 
Jp 

+8 fc(p) 2 v 2 H v (p)[B'(x v )-B'(xo)} 
Jp 

+ 16 / H v (p)[B'(x v ) ~ B'(xo)] 2 hip) 
Jp 

4 I [C(p)v 2 + Ahip) (B'(x v ) ~ B'(xo))} H v (p)[l + 2I 2 (p)H v (p)} 
Jp 

-32 fh{pf [C(p)v 2 + AI 3 (p) (B'(x v )-B'(x ))]H v (p) 2 B"(x ) 



+ 



N 



■16 / h(p)B'(x v ) (l + 2I 2 (p)B"(x ))H (p) 
Jp 

f \og[l- 4 I 2 (p)B"(x v )} 



+2 / I 2 (p)H v (p)[l + 2/ 2 (p)S"(xo)][l - 4/ 2 (p)5"(x 

J P 

+2 fl 2 (p) 2 H v (p) 2 [l-AI 2 (p)B"(x )} 2 
Jp 



H v {p) 



-16/ 2 J I 2 (p)B'(x v )[l + 6I 2 (p)H (p)]B'"(xo) 
B"(X V ) 



(6.2) 
(6.3) 



l-4I 2 (p)B»(x v )' 

where v 2 stands for v 2 b and is the argument of the function. The integrals I n and I n (q) are defined in 
appendix I and we recall that: 

X v = v 2 + 2Th + AI 2 {B\x v ) - B'(xo)) (6-4) 

6.2 The effective action in terms of the renormalized disorder 

To make the conversion to renormalized disorder, we use that on the r.h.s. of (6.2) we can express B in 
terms of the renormalized B at leading order. Using (3.8) 

B'(x) = B'(x + 2Th + 4I 2 (B'(x) - &(0))) (6.5) 
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x :— v 
and differentiating (6.5) w.r.t x, we obtain 

B"(x) 



l + AI 2 B"(x 
This allows to rewrite (6.3) as (already noted in (4.25)) 

H x {p) :-- 



B"(x + 2Th + 4I 2 (B'(x) - B'(0))) 



l + 4(I 2 -I 2 (p))B"(x) 



H v (p) 



(6.6) 



(6.7) 



(6.8) 



SB(x) 



N 



C{p) 2 H x {p)x 



+64 I 3 (p)C(p)H x (p) 2 x B'(x)-B'(0) 



+128 / h{pfH x {pf B'(x) - B'(0) 



C(p) 2 xH x (p) B'(x)-B'(0) 



+16 H x (p) B'(x)-B'(0) I 4 (p) 
Jp 



N 



C(p)x + 4J 3 (p) (B'(x) - fl'(O) ) H x (p)[l + 2I 2 (p)H x (p)} 



-32 / I 2 (pf \c(p)x + 4 J 3 (p) (b'(x) - B'(0))] H x (p) 2 D ' (0) 



1 + 4/ 2 5"(0) 



+ 



T 2 



16 / h(p)B'(x) I 1 + 2 Up) B "(°) 

l P K ' K '\ W l + 4/ 2 B"(0) 



1 f^(l + AB"(x)[I 2 -I 2 



Ho(p) 



N\ 2./ \ l + 4B"(x)I 2 



+2 / I 2 (p)H x (p) 



l + 2/ 2 (p) 



B"(0) 



+2 / I 2 (p) 2 H x (pf 



1 - 4/ 2 (p) 



l + 4/ 2J B"(0) 
B"(0) 



l-4/ 2 (p) 



B"(0) 



1 + 4I 2 B"(0) 



1 + 4/ 2 E"(0) 

B'"(0) 



-l^Jh(p)^)[i + ^(p)H^)] Ti 



(6.9) 



6.3 Expression in terms of rescaled dimensionless quantities 

As in Ref. [6], we define the dimensionless function b of the dimensionless argument z through 

b(z) = AA d m^- e B(zm- K ) . (6.10) 
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where e = 4 — d and A d = 2(47r)~ d / 2 r(3 — |). This prefactor has been chosen to simplify the leading-order 
/3-function. For d < 4 and T = all integrals which will appear here and below are convergent in the limit 
of UV cutoff A — ► oo that we consider, e.g. in that limit I 2 = A^m~ e /c. The rescaling factor ( is for now 
unspecified, but at the fixed point it will yield the roughness exponent. 

In terms of rescaled quantities, one finds that the result of the 1/iV expansion for the second cumulant 
of the disorder can be rewritten as 



b m (x)=b° m (x) + ^b 1 [b° m ](x) 



(6.11) 



where 6^(x) is the solution to dominant order, obtained previously for an arbitrary bare disorder. We have 
made apparent the dependence on the mass m (and raise the indices whenever necessary). The expression 
for the 0(1 /N) correction bi[b](x) can be obtained from the last subsection upon rescaling and reads (for 

d<4): 



bi[b](x) 



1 
A~ d 



+ 



A, 



2c(p) 2 [h x (p)x] 2 + 4t 3 (p)c(p)h x (p) 2 x [b'(x) - b'(0)} 

+2i 3 (p) 2 h x (p) 2 \b\x) - b'(0)] 2 + 2c(p) 2 xh x (p) \b\x) - 6'(0)] 
+ h x (p)[b'(x)-b'(0)] 2 U (p) 



-Mp)h x (p) 



[c(p)x + i 3 (p) (b'(x) - b'(0))] h x (p) 



2rp2 



+ 



4A" 



-U 2 (pf [c(p)x + i 3 (p) (b\x)-b\0))]h x (p) 2 ^^ 

^l ln fl + b"(x)[t 2 -t 2 (p)} 
2 n 



1 + b"{x)i 2 



1 



+-h(p)h x (p) 



+^2{pfh x {pf 
-i 2 i 2 (p)b' \x) 



1+ 2 12{P) 1+^(0) 



l-i 2 (p) 



&"(0) 



l + i 2 6"(0) 



«2(p) 



6"(0) 



l + i 2 6"(0) 



where one has defined 



hx(p) 



-i 2 (p)h (p) 



b"{x) 



6"'(0) 



;i + ^ 2 6"(o)) 3 



as well as the rescaled temperature: 



T 



1 + [l2 - 12(P)]6"(X) 

4A d m e , 



-T 



(6.12) 



(6.13) 



(6.14) 



where = d — 2 + 2£ is the energy fluctuation exponent, and the rescaled integrals are denoted by small 
letters: 

i r d d p 

m=i = A~ d J ^p 



*n(p) := -1-4 (p) 



(6.15) 
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c{p):=—- 2 (6.16) 

1 + p z 

with, e.g %2 = i 2 (p = 0) = 1/e. In the expression above we have kept the order of the diagrams from section 
4.3 ff.. Note that explicit A dependence can be reinstated in (6.12) by restricting all rescaled momentum 
integrals by K/m as upper cutoff, and is necessary in T > integrals (since they are usually UV divergent). 
Regrouping terms, this result can be rewritten in a more compact form: 

h(x) = 2x 2 gi (a x ) + 2x{b'{x) - b'(0))g 2 (a x ) + 2(b'(x) - b' (0)) 2 g 3 {a x ) (6.17) 

+eT m [x(g A (a x ) + a g 5 (a x )) + (b'(x) - b'(0))(g 6 (a x ) + a g 7 (a x )) 

+ b'(x)((g 8 (a ) + a g 9 (a ))] 
+ {eT m ) 2 [g w (a x ) + a Q gu(a x ) + alg 12 (a x ) + b'{x)a(^ + a g 13 (a ))] 



with 



b"{x) _ a x _ b'"(0) _ _ 6 ,/// (0)/e 26'"(0) 2 /e 2 

"x\P) , • / \ 5 Oi h"(0\\n ' ^ 



and 



i + ^ M , x^j i-a x i 2 (p)' (1 + 5^21)3' [l + 6"(0)/e] 2 [1 + 6"(0)/e] 3 

(6.18) 

9iM = ^~ fc( P ) 2 h x (p) 2 (6.19) 



A 



g 2 (a x ) = -^ f[2c(p) l3 (p)h x (p) 2 + c(p) 2 h x (p)} (6.20) 

A d J p 

g 3 (a x ) = -j- j[\k{p)h x ip) + i 3 (p) 2 h x (p) 2 ] . (6.21) 

The other functions gi(a) and 7, which characterize non-zero temperature, are given in Appendix H. 

Finally, note that the results (6.12) and (6.17) given above are for the choice ( = 0. For a non-zero 
C = Co + ^Ci + --- the result for b\ [b] (x) is identical to (6.12), (6.17) above, up to the trivial linear rescaling 
term 

61 [6] (x) -»• h [b] (x) - (2C 1 xb'(x) + 4Ci6(x)) ln(m) . (6.22) 

7 The /3-function at order 1/N 

We are now ready to obtain the flow equation of the dimensionless disorder, i.e. compute the /5-function. 
Here we give the most direct method to do so and give the results. A second method, closer in spirit to 
the diagrammatic approach and which yields more compact expressions has also been devised. Since it is 
rather involved it is detailed in appendix G.2. 

Our goal in the present paper is as follows. The dimensionless disorder b(x) depends on the IR cutoff 
m, and a priori, also on the UV cutoff A. To obtain a FRG flow equation we want to express: 

-md m b(x) = P[b](x) = (3o[b](x) + ^Pi[b](x) + ■■■ (7.1) 

in terms of b(x), at fixed A. Furthermore, we are interested in the behavior of the resulting expression when 
m/A becomes very small, which we hope can be made independent of K/m, if necessary with appropriate 
redefinitions of parameters. 
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To compute /3[b] (x) in the l/N expansion we write 

b(x) = b (x) + ^bM(x) + O (Jp) , (7.2) 

where b (x) is the dimensionless disorder at leading order. The corresponding /5-function was derived in 
Ref. [6]. It can be recovered by inserting (6.10) into (3.14) using md m I\ = —2m 2 l2 and the above value 
of I 2 (for A = oo). In the variables b it reads 

-md m b (x) = Po[bo](x) 

Po[b](x) = (e - 4C)6(x) + 2(xb'(x) + h'(x) 2 - b'(x)b'(0) + T m 6 ' ( ^ , (7.3) 

where we recall T m = AT A d m e / e. Thus the /3-function at leading order has a simple and well-defined 
A = oo limit. 

We now turn to the next order correction. 

7.1 The /^-function at T = 

We first detail the T = limit, which has a well-defined A = oo limit. At order l/N, deriving (7.2) w.r.t. 
m, we obtain 

- m d m b = p [b } + ^p [b ]^^ + o(^j . (7.4) 

Replacing b by b, using (7.2), we obtain the /3-function at order l/N 

= A> 



b-^hlb] 



♦^w¥+°GW 



This expression is still symbolic. The derivative ^6i[6] e.g. is the sum over the derivatives w.r.t. all deriva- 
tives of b, and in the above expression is multiplied by the /3-function of the corresponding derivative, 
obtained by deriving (7.3) w.r.t. x. Since at T = 0, bi[b](x) depends only on b'(x) — b'(0) and b"(x), (7.5) 
gives: 

« - s -ww M " ] " (x) + mw%) mb] ' (x) ~ mm 



e6i[6](x) + 6i[6]'(x) [ft'(x) - 6'(0)] - 6'(x)6i[6]'(0) [> (7.6) 

Inserting the expression (6.12) with T = one finds: 

/3(6(x)) = eb(x) + -b'(xf - b'(x)b'(0) ■ ' ' l <[ "'" 



2 w v / v / A^J (2vr) rf 

2a;c(p) 2 /i a; (p) 2 [a;e + 2xeh x (p)i 2 (p) + 21/(0) - 2b'{x)\ 

+4c(p)h x (p) 2 i 3 (p)[b'(x) - 6'(0)][2xe + 2xeh x (p)i 2 (p) + b'(0) - b'(x)} 
+2eh x (p) 2 [3 + 2h x (p)i 2 (p)]i 3 (p) 2 [b'(0) - b'(x)] 2 
+2c(p) 2 h x (p)[b'(x) - b'(0)][xe(l+i 2 (p)h x (p)) + 1/(0) - b'(x)\ 

+eh x (p)[2 + h x (p)i 2 (p)]h(p)W(0) ~ b'(x)} 2 } . (7.7) 
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A more compact expression can be found if one uses (6.17) as a starting point. The first term on the l.h.s. 
of (7.6) is replaced by J^-(— md^ax) and one uses 

-md° m a x = ea x + [b'(x) - b'(0)]b"'(x)(a x ) 2 = ea x + [b'(x) - b'(0)]d x a x . (7.8) 

One finds a form rather similar to (6.17): 

1 



-md m b = eb + ^b' 2 - 6'6'(0) 



+ N 



2x 2 g 1 (a) + 2x(b'(x) - b'(0))g 2 (a) + 2(b'(x) - b'(0)) 2 g 3 (a) 



(7.9) 



with 



9i( a ) = e (ag[(a) - g^a)) (7.10) 

g 2 (a) = eag' 2 (a) - 2g 1 (a) (7.11) 

93(a) = e(g 3 (a) + ag' 3 (a)) - g 2 (a) . (7.12) 

All terms proportional to d x a x have canceled. Also note that we have used analyticity in the derivation. 
Issues related to the non-analytic regime will be discussed in a subsequent publication. 

7.2 The /^-function at non-zero temperature 

At non-zero temperature the /3-function to order 1/N is independent of the UV cutoff A only for d < 2. 
Its expression is more complicated and we give here only its form, for d < 2. The derivation and explicit 
expressions for the functions are given in Appendix H, together with some comments about d > 2. 

j3 1 [b](x) = f3l =Q [b](x) + T m [x(g 4 (a x ) + a g 5 (a x )) + b'(x){g 6 (a x ) + a g 7 (a x ) + g 8 (a ))} 

+T^ [<7io(a x ) + ao<7n(a x ) + alg 12 (a x ) + xa' (e - a )g 5 (a x ) + ab'(x) {eg 7 (a x ) + <f>(a ))] 

+T? n b'(x)(ij(a )a 2 + 4){a )a) + ae 2 (g u (a x ) + 2a g 12 (a x )) (7.13) 

where the T = expression Pf =0 [b](x) was given above. 

8 Conclusion 

In this article we have computed the effective action of the field theory of random manifolds at large N. 
The 2-replica part of this quantity is what is needed to compute the renormalized disorder to order 1/N. 
Although similar in spirit (one must compute the determinant of fluctuations around a iV = 00 saddle point) 
the problem solved here is much more complex than for the standard 1/N expansion (say in the 4 -model), 
first because one needs to perform the calculation of fluctuations around the saddle point at fixed averaged 
field value, second because this involves four-replica matrices. It does however not involve spontaneous 
replica- symmetry breaking of the Parisi type, but rather some type of simpler explicit (vector) symmetry 
breaking. 

To handle such additional difficulties we have introduced in this paper two complementary methods. 
The first one is graphical and uses a diagrammatics which is able to handle both the O(N) and the replica 
indices. In this diagrammatics the zero-temperature diagrams are reasonably easy to compute. Much more 
subtle are the finite-T diagrams. Interestingly, only order T and T 2 are found to be non-vanishing to this 
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order in 1/iV. It is even simpler for d — 0, where the result for the /3-function is polynomial. (This can 
be derived by using Bogoliubov's subtraction operator.) The second method uses an algebraic formula for 
the determinant of fluctuations around the saddle point. The algebra of the four-replica matrices is worked 
out and one uses an expansion in number of free replica sums to compute all components iteratively in a 
given order of this expansion. Since we only need the two-replica part, this is a finite calculation, although 
rather tedious and has to be performed with Mathematica. The two methods are complementary, since 
any forgotten diagram of the graphical method can be traced to some term in the algebraic result, and vice 
versa. They are also, to our knowledge, new, henceforth the detailed exposition. 

Having obtained the effective action, we rewrote it in terms of the dimensionless renormalized disorder. 
By varying with respect to the infrared cutoff, we obtained the /3-function to order 1/iV. We noted that this 
/3-function is UV finite at T = 0. 

It is important to note that the derivation was made, strictly speaking, using an analytic action. This is 
familiar for N = oo, where the same strategy was applied successfully: Although the derivation was done 
in the analytic regime, the /3-function could then be continued to the non-analytic one. This was done via a 
careful analysis of the solution when it reaches the Larkin scale. A similar analysis will be performed in a 
forthcoming publication, together with a comparison to the two-loop result, and a detailed analysis of the 
physical consequences of the FRG flow derived here. 
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A The large-iV formalism for the effective action 

Before studying specific models, in the following sections, we first present schematically the framework of 
the large- N calculations. 

A.1 General properties of the l/N expansion 

The general problem can be formulated as follows. We want to compute the effective action T[u] defined 
as the Legendre transform of W[J] = \nZ[J] (J being the source field conjugated to u), in the case where 
the partition function can be written as: 

Z[J]= [ V[ij}e- NS[ ^ j] (A.l) 



where j = J/y/N is the rescaled source and ip is some auxiliary field (or a set of such fields). Here 
all space coordinates and indices are suppressed and integrals and sums implicit, in order to exhibit the 
structure more clearly. 

The first step is to write yV[J] in an 1/JV expansion using the standard saddle point method. One finds: 

W[J] = NW[j] (A.2) 

W[J) = W°[j] + ^W l \j] + -^W 2 [j] + . . . , (A.3) 
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where the j-dependent saddle-point value ipj of the auxiliary field is solution of 

sifei] = ' ( A - 4 ) 

The expansion yields 

W°\j] = -S[rJ> j ,j] (A.5) 

W 1 [j} = -hr\n S';^,j] (A.6) 

W 2 [j] = —SZMjJ] {faMJt)^ + ^SZ c [^,j}S7 fg [^,j} (^aMcMM^-i ■ (A.7) 
More generally, the W n (j) are obtained from the loop expansion of the field theory: 

\nZ\j] = NW°[j} +lnj V\<j>] exp I --<t>S"\^ j]<f> - £ -^^ . ..d^StyjJ]^...^ 

' (A.8) 

In these formula, the indices a,b, ... summarize all spatial coordinates, indices etc. of the field. 
The effective action is then defined as: 

I>] =NT[v = u/VN} (A.9) 

F[v]=vj v -W\j v ], v = W'[ ]v ] (A.10) 

The equation for j v can formally be inverted into an expansion in 1/iV using (A.3), 

jv=fi + jjjl + --- , (ATI) 

which yields in turn the expansion for the effective action: 

oo 

f[ v ] =J2N- p T p [v} . (A. 12) 

One finds that the leading order is simply the Legendre transform of W^°[j], 

f>] = vj° v - W°[j° v } , v = (W°y\j° v } , (A.13) 

where here and below (W n )'[j] = ^"[j], Since one has c^T ^] = j®, it satisfies the self-consistent 
equation 

f°[v] = vd v F°[v] - W°(d v t°[v]) , (A.14) 



and one has the usual relation between the second derivative matrices (W°)"[7°] 
In this paper we use the result for the next order: 



6T° [v] 
SvSv 



T 1 [v] = -W 1 [£] (A.15) 
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For completeness we give the two next orders: 

f 2 M = -w 2 [f v ] + \{w')U\{d%wy\ji\-\w')U\ 

I5t l ( 8 2 Y G \ l 8V l 



^ + 2 *7 ( AA J AT (A ' 16) 



f>] = -w 3 ia + (w 1 )'Kl(»" > )"K]- 1 (w 2 )'b- 



01 



+l(wyy v ](wy\j° v }- 3 (w i y\j° v } 3 ca.iv) 

where the following graphical rules allow to restore correctly all index contractions and spatial integrals 
implicit in the schematic notation above. Denote — T n by a square with n inside and W n by a circle with 
a n inside. One treats the circled 1,2,3 etc as vertices and considers all tree graphs. A line is a propagator 
Wo '\jv\ 1 an d is thus not a vertex (it is summed in a line). The sum of these numbers is just the order. The 
formula are equivalent to: 



IdKD 




®-@ 4 (iHIHD 



This has, including the combinatorial factors (from expanding the exponential and number of choices), an 
immediate interpretation as a perturbative expansion with all 1PI graphs subtracted. In the last line, the 1-2 
graph comes with 1/2! in expanding the exponential, and adding the two possibilities, the 1 — 1 — 1 comes 
with 1/3 ! , but then there is a 3 for selecting the middle one, and the last graph has a 1/3 ! from expanding the 
exponential. There is also a relative minus sign from each vertex added. To justify this graphical method, 
one simply recalls that quite generally correlation functions, i.e. W[J] are given as the sum of possible tree 
diagram made with T[u] vertices. Since the Legendre transform is involutive, the same is true for T[u] in 
terms of Y\?\J] vertices. Thus to write T n [v ] one must simply insert the proper orders in 1/N at each vertex, 
in all possible ways, so as to match the total order. 

Thus all the T p [v } can all be expressed as functions of the W p (j®). Inserting the results (A. 5), (A. 6), 
a.s.o. from the saddle point expansion, one obtains the T p [v] explicitly in terms of derivatives of the S 
functional in (A.l). 

A.2 Summary of main result 

Before we detail the calculations in a more pedagogical way for some specific models below, we first 
summarize here in compact notations the main result for the two lowest orders of the 1/N expansion, with 
a generalization to a bilocal bare action. 

Let us consider a iV-component field theory whose action functional can be written as: 

S(<f>) = \(t>G- l ct> + NS- mt [^ xy \ (A.18) 
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V>: 



xy 



N 



'x ry ? 



(A.19) 



where S- mt is a functional of the bilocal field if) xy (which is also a bi-index matrix, if the field carries other 
indices (e.g. a,b, . . .). Then its effective action can be written as: 



1 



( t>G- 1 <f ) + NT ^ xy }+T 1 ^ xy } + . 



where T u is also a functional of a bilocal field and satisfies the self-consistent equation: 

5T° 5S- 

[ipxy] = Y7^ [V'xj/ + G [i/j] xy ] 



H. 



zt 



6if>. 



zt 



Gm 



xy 



G- l + 2 



5T° 



■ni 



and 



r 1 ^ 



■m\ 



-trln 



lzt,z't' + 2 



8 2 S, 



in I, 



■[ipxy + G[ij] xy ] X 



G[ip] ut 'G[ip] vz/ + iput'G[ip] vz > + i>vz'G[ip] ut < 



Where l xya hztcd — ^xz^yt^ac^bd- 



B Toy model 

For pedagogical clarity, we will give all details for the simpler case of the toy model. 

B.1 Model and effective action to leading order 

We study the following O(N) toy model, defined by the partition sum: 

Z[J]= fv[u]e- s M +VNju = fv[u]e~ s ^= f V[u\ V[ x \ V[X] e^"'*'^ 

= fv[ x ]v[\] e - NS[xXj] 

S[u,j] = J l -{Vuf + l -mW + NV ftp - VNju 



S[u,xA,J] = S(u) 



\Wx-u 2 ) 



5[ X ,A,j] = itrln(-V 2 + m 2 + A)+ I 



Nju 



V(X) - 2 A X 



(A.20) 

(A.21) 
(A.22) 



(A.23) 



(B.l) 
(B.2) 
(B.3) 



\ f J.(-V 2 + m 2 + A)^, (B.4) 

z J xy 



where X(x) and xi x ) are local fields (the factor i has been absorbed in the field A). The above expression 
(B.l) is thus of the form (A.l) where ip = (x, A) is a set of two auxiliary fields. 
The saddle point equation (A.4) thus read: 

Xj(x) = (-V 2 + m 2 + A,)^ + Jy ,{-V 2 + m 2 + A,)^(-V 2 + m 2 + X^jy (B.5) 

X j (x) = 2V'( Xj (x)). (B.6) 
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From Appendix A one finds the dominant order: 

W°[j] = -S[ Xj Aj,j] (B.7) 

Using the saddle point equation one finds: 

(W°)'\j] = -djS[ Xj , XjJ] = (-V 2 + m 2 + A,)" 1 ; (B.8) 

Thus one obtains: 

T(u) = NT(u/Vn) , v = u/VN (B.9) 

j v = (-V 2 + m 2 + \)v (B.10) 
f[v] =j v v + S\x Jv ,\j v ,Jv] 

= I,;(_v 2 + m 2 + \)v + itrln(-V 2 + m 2 + \ jv ) + NV( Xjv ) - \\ jv Xj v • (B-ll) 
One finds: 

f (v) = itrlnG; 1 + J l -{yvf + l -m 2 v 2 - l -(V 2 G v ) xx + V(v 2 x + (G v ) xx ) (B.12) 

(G- 1 )^ = [-V 2 + m 2 + 2V>* + {G V ) XX )]8 XX> . (B.13) 

This can be simplified further as the complicated part of the effective action is just a constant. First one 
notes that T(v ) = T[G V , v] satisfies: 

d G f[G,v]\ G=Gv =0. (B.14) 

As a consequence one sees that: 

sr\v] 



Svl 



(_V 2 + m 2 )v x + 2v x V (v 2 x + {G v ) xx ) (B.15) 

-'x 

The effective action (per unit volume Vt) for a uniform configuration of the field v x — v reads: 

If <„) = in.V + W + G.) + i ^Intf + m* + 2„V + G .)) + - - -^ ^ 2 - g - (B.16) 

G " = J, qi + mi + 2V V + G„) ' (B ' 17) 

which can also be written as 

Here one can write: 

f(v)=f[G v ,v 2 ] (B.19) 

d G f[G,v 2 ]\ G=Gv =0. (B.20) 

For uniform configurations the effective action is simply a function of v 2 such that 

\m 2 + V\v 2 ) = 1^1 = im 2 + W + G„) (B.21) 

i.e. it satisfies the self-consistent equation 



-f (u) = imV + V> 2 + G„) - G v V'(v 2 + G„) + \ I ln(g 2 + m 2 + 2V'(v 2 + G v )) . (B.18) 



J (? g 2 + 2\/'(t; 2 )/ 
This defines a renormalized potential V(v 2 ), whose RG flow is studied in Appendix I of [8]. 

31 



B.2 l/N corrections to the effective action 

To next order one has, from (A. 6) : 



W 1 \j] = ~trhiS"\x J ,X j J] (B.23) 



(S'^ x )xy = V"(x(x))S xy , (S'^ x ) xy = --5 xy (B.24) 



where Xj> <\,- are taken at their saddle point values (B.6). The symmetric matrix of second derivatives reads: 

1 

— < 
2 

\S\\)xy — —-^GxyGyx ~ {jG) x G xy {Gj) y (B.25) 

G xy = (-V 2 + m 2 + A)" 1 . (B.26) 



It can be put in the form: 



S'^-Uil). (B.27) 



Its inverse reads 



(ST 1 = "2 n Ti W n /.W (B- 28 ) 



-5(1 -AB)" 1 (1-fiA)- 1 
l-AB)" 1 -A(l-fiA)- 

A XJ/ = -2y"(x(x))<J xw (B.29) 

■L*xy ^xy^yx T" ^\Jv^ Jx^xyy^J )y ^xy^yx ~r ^Vx^xyVy > ^d.jUJ 

where the last equality holds only for j = j°. Note that A and £> do not commute (unless fields are 
uniform). 

The result for the effective action is thus, from (A. 15) 

r 1 (v) = ±tr\nS"( Xj °A j °J° v ) (B.31) 

3° = G~ l v (B.32) 

G xy = (-V 2 + rn 2 + \ 3 -o)£ (B.33) 
A more explicit form can be given for a uniform field configuration v x = v, with in that case: 

A xy = -2V"(x)$xy , X = Xj° = v 2 + G xx (B.34) 

B xy = (G x /) 2 + 2v 2 G? (B.35) 

,/,W^_ 2 fU v (q) + 2v 2 G(q) -1 



{S " ] iq) l + 2V"(v" + J k G v (k))[U v (q) + ^ 2 G(q)]{ ' "1 -2V"(v 2 + J k G v (k)) 

(B.36) 



n,(g)= / G v (q-k)G v (k) (B.37) 

^ = jfe2 + m 2 + 2 y^2 + J fcGi;(jfc)) ( B38 ) 



yielding finally 



r\v) = \ Jin (1 - i4(g)B( g )) = ^ J In (l + 2V" (v 2 + J G v (k)\ (U v (q) + 2v 2 G q {v))\ (B.39) 



up to a constant 
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C Calculation of order l/JV for the random manifold 

We now sketch the derivation of T 1 for the case of the random manifold. 

It was shown in Ref. [8] that the partition sum (2.4) can be put in the form: 

Z[J]= f D X D\e~ NS[x ' x ' j] (C.l) 

S[ X , X,j] = itrln^- 1 + iX) + J U( X (x)) - l -X a \x) X a \x) 

~\ I JaWiC- 1 + iX)- a l bx , 3b (x') , (C.2) 

^ Jx,x' 

where the inversion and trace are performed in both replica space and spatial coordinate space. It has again 
the form (A.l) where ip = (x ab ( x )i X ab (x)) is a set of 2-replica-matrix auxiliary fields. The saddle-point 
equation (A. 4) reads (see Ref. [8]): 

X f{x) = (G 3 ) axM + {G j : j) ax ■ (Gj : j) bx (C.3) 

tXf(x) = 2d ab U( Xj (x)) (C.4) 

Gj 1 = C- 1 + iXj , (C.5) 

where Gj is a matrix with both replica indices and spatial coordinates and inversion is carried out for both. 
Here and below, replica indices are raised whenever explicit dependency is given, e.g. Xab = xT- The 
notation for the iV-component vector (G : j)\ x = ^ c J G bXtCy Jl(y) is a shorthand for a matrix product, 
and everywhere we denote by 

d ab U(cf>) := d^JJ{<j>) (C.6) 

the simple derivative of the function U(4>) with respect to its matrix argument 4> ab . 

The dominant order, W°[j] and T°[v], was computed in Ref. [8], thus here we study only the next order. 
It is given by: 

T 1 [v] = Ur\nS"[xf,Xf,j } (C.l) 

fax = (G" 1 : v) ax (C.8) 

(G-^a^y = (C-^ySab + 2d ab U( Xv (x))S d (x - y) . (C.9) 

Note that when computing the fluctuations around the saddle point we consider Xab and xta as indepen- 
dently fluctuating fields, symmetry being restored at the saddle point (it is also possible to perform the 
calculation with symmetric matrices only). 

The matrix of second derivatives can again be put in the form (see the previous section): 

and one can show that: 

trlnS" = trln(l- AB) (C.ll) 

A ab x,cdy = -2d Xab d Xcd U(x(x))8 xy (C.12) 

^abx,cdy *-* ax ,dy^ cy ,bx > \J • *-* Jcy^aXjdyxS 1 • J )bx "T \J ■ V-T jax^bx^yy^ • J )dy t ^J-ab,cd 

L 7 ax ,dy*~J cy ,bx "T V bx \j ax,dy^cy \ ^ ax^J bx ,cyV dy ~r ±vlab,cd \\..lj) 
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and note that A and B do not commute (unless fields are uniform). The last equality is valid only when 
j = j® is inserted. 

For uniform fields one finds: 

f 1 [v] = ^jMl-A(v)B q (v)) (C.14) 

A(v) ab , cd = -2d Xab d Xcd U( Xv ) (C.15) 

B q (v) ab , cd = f G a v \k)G cb (q -k) + v b v c G a v \q) + v a v d G b v c (q) (C.16) 

Jk 

Gl\k) = [C(ky l +i\ v ]- b l (C.17) 

Xl d = v c v d + JGf{k) (C.18) 

Jk 

i\ ab = 2d ab U{ Xv ). (C.19) 

We note that the above matrices are not a priori symmetric in a <-> b or in c <-> d since we have chosen the 
representation where all components of the fields fluctuate freely. The alert reader will notice later (when 
coming back to the main text) that in fact the final matrix AB q will possess such a symmetry. 

D The excluded replica-formalism 

In this section, we calculate the terms at order T and T 2 using the excluded replica-formalism. This gives 
an independent derivation of the moon-shaped building blocks ^=^ and /=^ in section 4.2. 

We start by recalling that corrections at finite temperature are more difficult to obtain. The simplest 
example is given in the next appendix E. Here we want to understand this by making two contractions 
between B(u a (x) —u b (x)) and B(u a (y)—u b (y)). Be the first contraction (focusing on the replica- structure) 



■ ■ . (D.i) 

The following possible contractions are 

. . ; ; " i--"^ . " . ^~~~~~-; • -' 

At zero temperature, we have found in section 4.3 that only the first term contributes. This enabled us to 
do the calculation at T = 0. 

The general case is more complicated, since there are now three more possible contractions, one may 
draw, i.e. 4 for each B" . This looks discouraging. Two ways out of the dilemma can be thought of. The 
first, called recursive construction or successive construction, tries to add one more link to a chain of B". 
The decisive simplification then is that whenever arriving at e.g. 

(D.3) 

no more links can be added to the right, since the rightmost vertex is B" a , which is a constant. This 
procedure has however one crucial deficiency: Contractions are non-commutative! Let us illustrate this on 
the diagram 

If one first does the two leftmost contractions, then the two rightmost contractions are not possible. On 
the other hand, when doing first the two rightmost contractions, the two leftmost are possible! Working 
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with this formalism, i.e. using its implied simplifications, one has to number the lines. This is pretty 
awkward, since one wants to be able to interpret diagrams as such, without having to number lines. This 
is especially important here, since one wants to recognize the chains which give rise to H a b(p). Also note 
that problems arise at finite temperature, since there one has additional lines to add, and not all lines are 
needed to get a 2-replica term, i.e. there are lines which can be discarded, and still the diagram would be 
a 2-replica term. A very instructive example is the simplest 2-loop order diagram at finite temperature, 
which is derived in appendix E. One sees that diagrams can be grouped differently to cancel, and that these 
different cancellations correspond to different paths in the contractions. 

One way out of the dilemma is the excluded replica approach, which at 2-loop order has as its descen- 
dent the sloop-algorithm [4]. Here we explain the excluded replica approach. We start by writing the trivial 
relation 

J2 Bab = J2 Bab + Y, Baa - ( D - 5) 

a,b aj^b a 

The first term on the r.h.s. is the excluded replica disorder, the second a constant. When performing contrac- 
tions, the second term does not contribute. Thus one can perform perturbation theory instead with J2 ab B a b 
with XL^b B a b- The big advantage is that the last two terms in (D.2) do not contribute. The backdraw is 
that the final result is a sum over excluded replicas, which has to be projected onto the 2-replica term, and 
one may have more terms at intermediate steps. The final projection can be done formally, by replacing 
B a b — ► B a b(l — Sab), expanding and then collecting the 2-replica contributions. 
We now introduce the excluded replica-formalism. Recalling that 

JU I O I = i _ 4/ 2 ( p ) B"( X ab) * ^ * = 2 ^ Ua ~ Ub)21{ HM ^ Ua ~ ^ ' (D ' 6) 

We then proceed as follows: Whenever we have a long chain with some replica conserving lines (as above) 
and with doubled lines (or single lines when inserting vertices), we will always sum the double-line part, 
resulting in H ab . We then have to sum explicitly the rest, and project it onto the 2-replica-contribution (in 
general). For illustration, we start with the open chain. This is, say when fixing one end to have replicas a 

and b 

p a 



We furthermore observe that the combinatorial weight for either having a replica-conserving double line 
( <^> ), or a doubled line up ( ■'O* ) or down ( «cz>" )■> * s me same. We do not draw crossed lines; 
by definition they are incorporated into the combinatorial factor. The above has to be projected onto the 
2-replica sum. Writing any H c d{p) with excluded indices as H c d(p)(l — S c d), and multiplying all terms, 
one always has to take the term with the 5 c d, and (D.7), summed over all lengths of the chain, is 

oo 

E -^ l OoOo<yxxx> = H ab (p) x>4/ 2 (p)# aa r 

all lengths «=0 

= H ab (p) l 



l + 4I 2 (p)H aa (p) 

= H ah {p) (1 - U 2 {p)Bl) . (D.8) 

This object is well known: It is the chain times the moon-diagram, see table 5. We remark the important 
point that, even though we have an infinity of diagrams to sum up, in the projection onto 2 replicas (or since 
we are already fixing two replicas, maybe we should better say onto one replica) all but 1 non-trivial term 
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vanish. (Also remark that we have been a bit sloppy to amputate the last u 2 legs. Otherwise, that is not so 
easy to represent.) 

We need another intermediate result. Summing all diagrams which conserve the lower index, we obtain 



We want this sum 



= <z> +...+ 
projected onto 1 replica, 



(D.9) 



Using the same procedure as above, we find 



S(p) 



oo 



2J 2 (p)2J-2J 2 (p)if a 

n=0 

2/ 2 (p) 
2/ 2 (p) 



1 + 2I 2 (p)H aa 
l-4/ 2 (p)i?l 
1 - 2I 2 {p)B'l a ■ 



The 2-replica contribution is 



o 2 



x 



S(p)H ab (p)S(p) . 



(D.10) 



(D.ll) 



This half-chain is very practical, since when inserted into a more complicated diagram, no further restric- 
tions apply. In a diagram involving this half-chain and some "rest", there are two 2-replica contributions: 
The half-chain projected onto the 1 -replica part times the rest projected onto 2 replicas (taking care of the 
replica-conserving line of the half-chain) or the 2-replica part of the half-chain times the 1 -replica-part of 
the rest. 

We continue on the order T 2 -term. We have the following terms before projection: 




... . (D.12) 



A dashed line denotes the identity; one would like to print this on a torus. The first term is very symbolic, 
since it is not a geometric sum but a log, and there should be at least one vertex. The second is (coinciding 
replicas in the chain) 

^ (D.13) 



0. 



We now have the choice of where to project onto the 2-replica term. The rest has to be projected onto the 
1-replica term. We first project one of the half-chains onto a 2-replica-term, using (D.ll). The 2-replica 
term is shaded in dark grey. The rest consequently has to be projected onto the 1-replica terms, shaded in 
light grey. We obtain the class C\ : 



Ci 




(D.14) 



Denote by n the number of H ab in each diagram, which have lines starting at different replicas. Each such 
H ab contributes a factor of (1 — Sab), for a product of 



E 



(1 - S hi2 )(l - S i2i3 ) ... (1 - S inh ) 



(D.15) 
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We need all indices to be restricted to be identical. The terms which lead to that restriction are either the 
product of all <Ts (1 term) or any one of the 5's left out. This gives another n terms, but with a different 
sign, for a total of 

OO 

& = S(p)H ab {p) J2(l ~ n) [-H aa (p)S(p)} n , (D.16) 

n=0 

where it is important to note that also the first term in (D.14) for n = is correctly given by the above 
formula. We note the auxiliary sum 

We give the intermediate results 

\-2I 2 {p)B" aa (D.18) 

21 2 {p){\ - AI 2 {p)Bl) (D.19) 



1 + H aa {p)S{p) 
S{p) 



the above sum is 



l + H aa (p)S(p) 

i+2 ^^- l-Sz (d - 2o) 



Ci = 2I 2 (p)H ab (p)(l + 2I 2 (p)Bl)(l - 4I 2 (p)Bl) . (D.21) 



This reproduces diagram (4.52). 

We now turn to the diagrams where the half-chain is always projected onto the 1 -replica contribution, 
whereas the non-trivial terms come from omitting a 5 belonging to one of the H ab with lines entering into 
different replicas. With the same shading for 1- and 2-replica terms as above, this is 



C 2 = (^ \ + ,,"^ V^vz^n + ,,'^T V^\ztV^- \ t +... (D.22) 
The first term is 0, since it necessarily has crossed indices: 

i =0. (D.23) 



Each H a h which is projected onto a single replica again contributes a factor of (1 — 5 ab ) Starting from the 
second diagram, we have to leave out exactly 2 8's in 

±(l-6 ili2 )(l-5 i2i3 )...(l-5 inh ). (D.24) 

n 

The first has already been left out in plotting (D.22), accounting for the factor of n. Leaving out two 8's 
leads to a factor of n(n — l)/2, since they are indistinguishable. The result contains two factors of H ab ; the 
remaining factors are all H aa : 

1 OO 

C 2 = -H ab {pfS{pf J> - l)[-H aa (p)S(p)T- 2 . (D.25) 

n=2 
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With the sum 



JT(n-l)(-x) 



n-1 



1 



n=2 



(1 + x) 2 ' 
and using (D.10) and (D.18) we obtain 

C 2 = 2l 2 {pfE ah {pf{\-U 2 {p)B" aa f 



(D.26) 



(D.27) 



This reproduces the term (4.53). 

We now turn to the diagrams of order T. First of all, note that the following two diagrams have no 
contributions proportional to B'^ a : 



and 




(D.28) 



The reason is that they have no doubled line, as in 




and 




(D.29) 



Note that both diagrams necessitate a replica-conserving double line (i.e. »<3>» or «<3>» ). Otherwise 
they vanish. The replica-conserving double line can be more generally replaced by all chains, which start 
and end with a double line, and which conserve the index running through. The chains with one index 
entering and one index exiting are 



Each of these has the form (with n + 1 half-chains) 

n 

S(p)H((l-6 ijij+1 )S(p)H ijij+l (p)) 



(D.30) 



(D.31) 



i=i 



The only replica-conserving term is obtained by using 5j.j. +1 in each factor. The projection onto the 1- 
replica term is denoted 1Z (for replica-conserving) 



K = 
and is evaluated as 



S{p) 



n = Sip) J2 l-S( P )H aa (p)] n = = - q( { Z . ; = 2/ 2 (p) (1 - AI 2 {p)Bl 
t^ l+S{p)H aa {p) 



... (D.32) 



(D.33) 



?1=0 



where the last identity can be found in (D.19). This reestablishes the two factors of (1 — 4I 2 B" a ) in diagram 
(D.29). Note that 1Z is /^=^ introduced in section 4.2. 

We can also give an equivalent interpretation of (D.21). It is H a b times summed half-chains, but since 
the indices are forced to be equal at the end, we can drop one of the S's in (D.3 1). Instead of (D.33) this is 



T = S<p)J2[-S<p)H aa <p)] n (l 



n 



n=0 



2I 2 (p) (1 + 2I 2 {p)Bl) (1 - 4I 2 (p)Bl) , 



(D.34) 



where the combinatorial factor (1 — n) is due to the fact, that one can drop one of the S's. We have also 
used the simplifications of equations (D.17) ff. T is nothing but /=x , introduced in section 4.2. 
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E The 2-loop diagram with a tadpole and graphical interpretation of 

perturbation theory 



For simplicity of notations, this calculation is done for a 1-component field u. We also note R(u) := B{u 2 ). 
We here calculate the 2-loop diagrams at finite temperature. This shows how the naive rules one uses at 
zero temperature can be misinterpreted. 

All diagrams have two vertices R, two lines between these two R's, and a tadpole attached to one of 
the i?'s: 



2- 9-. 9-, 9- v-, 9-, 9- 9-, 9-, 9- 



.-9 J. 



J/ 



-9 



-9 .-9 <? 



. i 



(E.l) 



This is agraphical representation of R'^R'^ d (S aa +Sbb — 28 c u ) ) {&achd + ^adhc) , together with combinatorial 
factors. Projecting onto 2-replica-terms only gives: 



9-, 



?~. 9-, 9-. 9- f 



.-9 



-9 .-9 _9 ... 



(E.2) 
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The interesting diagrams proportional to R""(0) are: 





(E.3) 

One sees that one has cancellation either in one column or in one line. One can interpret this as follows: 
Making first contractions between the two i?'s, identifying the replicas of the upper R, one can say that no 
tadpole can be added to the upper R. This is the cancellation in the first column. Equivalently, one can first 
draw the tadpole. If one does this connecting replicas a and b on the upper R, then one can say that it is 
no longer possible to add correlators connecting the two i?'s. This is the cancellation in the first line, as 
opposed to the first column. 



F Details of the algebraic method 

F.1 Algebra of 4-replica matrices 

The 4-replica matrices needed in the tr In calculation of Section 5 can be parameterized as in (5.5) and 
(5.10). They form an algebra, and to solve Eq. (5.11) one needs to write explicitly the components of an 
arbitrary product of such matrices. 

Let us consider two matrices M 1 and M 2 parameterized respectively by (x a , y ab , z ab) t abc) u l ab , v abc , w l ahc , 
g abcd ) for z = 1,2. Then the product M ab:Cd = Y, e f M ab,ef M ef,cd is itself parameterized by (x a , y a b, z ab , t abc , 
Uab, Vabc, w abc , g abcd ) and one finds: 



2 
afd 



X a — X a X a + 2y aa X a + 2x a Z aa + 2 2_^ Vaf Z af 

f 

Vad = x a y 2 ad + 2x l a t 2 aad + 2yl a y 2 ad + 2 J^ y\ f i 

f 

z ab = z ab x a + 2t aba x a + 2z ab z aa + 2 2 ,t abf z af 

f 

t a bd = z ab y ad + 2t aba y ad + 2z ab t aad + 2 2_^ t ab ft a f d 

f 

' X a U ac ' U ac X c ' V acc X c ' X a W aac ' ^Vac Z ca ' ^Vaa U ac ' ^ U ac Z cc 



(El) 



u 



^ V aefWefc 



+ Yl [ 2 vlf W lfc + 2v acf z lf + <fU% + vl ff U% + ul f W 2 ffc ] + 

f ef 

v ac d = x a v acd + x a g aacd + 2y ac t cad + 2y ad t dac + u ac y cd + u ad y dc + 2y aa v acd + v acc y cd + v add y dc 

+2 U a C t C cd + ^ U ad^ddc 
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Er -| n 10 1 O 1 O 10 I 1 ^ ^ 1 

[^Vafdafcd + 2 V acft c fd + 2v adf^dfc + U af V fcd + U af9ffcd + V aff V fcd\ + ^, V aef9efcd 
f ef 

Wabc = W abc X c + 9abcc X c + ^abc Z ca + ^bac Z cb + Z ab U ac + Z ba U bc + Z ab W aac + Z ba W bbc + 2"W afec Z cc 



,2 



+2£ a6a 'U ac + 2t bab ul c 



+ Yl [ 2 9abcf Z cf 2t lbfW 2 afc + ZtbafWbfc + W lbfU 2 f c + #a6// M /c + W afe/^//c] + ^ 9lbef™lfc 

f ef 

_ 1 2,1 2,1 2,1 2,12 ,12,12 i 12 

SUcd — W abcVcd ' W abdVdc ' 9 a bccy c d ' 9abddVdc ' Z ab V acd ' Z ba V bcd ' Z ab9aacd ' Z ba9bbcd 

^^abJ-cad + ^baJ-cbd + 2 ^abdPdac + ^bad^dbc + ^abaPacd + 2 Hab V bcd + 2w ab<$>ccd + 2w abdPddc 
+ ^, [ 2 ^abf9afcd + 2 ^baf9bfcd + 2 9 a bcf^cfd + 2 9abdf^dfc + W abf V fcd + 9abff V fcd + W abf9ffcd\ 



E 



9abef9 e fcd 



where we have made replica sums explicit. 

Using these multiplication rules one can rewrite Eq. (5.11) in terms of a set of nonlinear equations 
for the components of M A in terms of the components of .M. Unfortunately no closed solution seemed 
possible (except in some very special cases). The next step is thus to expand each component in number of 
replica sums as in (5.13). This results in a hierarchy of equations for components with increasing number 
of replica sum. For instance the zero-sum components behave under multiplication as in (F.l), dropping 
all terms with replica sums, and so on. These equations can be solved iteratively, as discussed in the text. 
For this, one needs the zero- and one-sum components of the matrix M., the calculation of which we now 
detail. 

F.2 Calculation of the matrix M 

We start by computing the matrix M in (5.1). As in the following subsections there are two stages. First 
make all Kronecker-deltas explicit, then expand each term in the number of replica sums. As discussed in 
the text, all intermediate free sums over two or more replicas can be dropped. 

A straightforward calculation from (2.6), for a model with only a bare second cumulant, gives: 



M abtCd = (2Td Xab d Xcd U(x))\ x = Xv (F.2) 

2 



-—{Sabcd 2_^ B" e + b ab b cd B" ac — (5 acd + 5 bcd )B'^ b — (5 abc + 5 abd )B'{. d + (5 ac 5 bd + 5 bc 5 ad )B" b ) 



where B" ah = B"(xf)> and x ab = xT + xl b ~ 2x" 6 ; Xv b * s gi ven by its saddle point value (2.16). Note 
that since one takes everything at the saddle point at the end, which is symmetric in a, b, all expressions 
resulting from two derivatives are symmetric: B" b = B ba (even if one chooses the fluctuating fields a priori 
non- symmetric). 

The matrix x% b still contains explicit Kronecker-deltas. As in the main text, one writes 

X? = Xab + SabXa , (F.3) 

where x ab and Xa contain no Kronecker delta, and are computed below. Then one sees that 

B"(X?) = Sab [B"(Xaa + Xa) ~ B" (Xaa)] + B" (^ab) • (F.4) 
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Inserting this form into (F.2) above one finds that the contribution of the piece 5 a b[B"(xaa + Xa) — B"(x aa )] 
cancels exactly and thus one obtains that M is given by (F.2) but with now B" ab = B"(x ab ), i.e. the part 
with no Kronecker delta. 

We can now continue the calculation from (2.16) by expanding in the number of replica sums. First we 
define: 



-2T3U 



Sab 2J U a 



u, 



ab 



u nh = -Lb' 



T 






(F.5) 
(F.6) 



where B' ab = B'(vl b ) and similarly for the three-replica term (which will drop later on). From there and 
(2.16) we obtain, dropping all higher order sums, the expansions: 



Xa 



■2Th-2ThY.Ua 



Xab = v 2 ab + 2Th + 2TI 2 



U ab - AUaa + U bb ) + - J](C/ ae + U be ) 



(F.7) 

(F.8) 



+2T/ 3 



J2 U ab{Uae + U be ) - J^ U aeU e b ~ J2^ UaaUae + UbbUbe _ o^^ea ~ -U be U eb ) 



;-» — LLC ~ CU (-» 

Thus in each of the B" ah = B"{x ab ) matrices appearing in (F.2) the argument can be Taylor expanded, i.e. 

B: b =Xb + BlJ2 ^fi CF.9) 

/ 



as 



where we have defined B ab = B"(x a b), B ab = B m (x a b) and 



1,*. 



Xab = Vab + 2T/x + 4I 2 (B' ab - -(B' aa + B'J) 



4 1 ~, 

Oabf = j;h(S' abf - -(S' aaf + S' bbf )) + 2I 2 (B af + B' b 



2 
bbfj 



bfJ 



+7fh 



B'ab(B af + B' bf ) - B' af B' ef - (B' aa B' af + B bb B' bf - -B' af B' fa - -B' bf B' fb ) 



(F.10) 



(F.ll) 



It will turn out below that at the end we will only need O aa f = AI 2 B' a f. We will not perform this expansion 
and replacement now. First, we turn to the calculation of the matrix N and perform the product MN , 
keeping B" ab unspecified. 

F.3 Calculation of the matrix W 

We now compute the second matrix, N , expanded up to one free replica sum. One has: 



N lb,cd = v a v d G q bc + v b v c G q ad + Tlf abcd 
G q = C(q)5 + J2 C(q) n+1 (-2TdU°) 



(F.12) 
(F.13) 



n>l 
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Kb,cd = J{Ad5 bc + J2 Jin + i(Sad(-2TdU°) bc + 5 bc (-2TdU°) ad ) 

n>\ 

+ E jq rn + i,n + i(-2TdU°)Z(-2TdU°)l (F.14) 

m>l,n>l 
f 1 1 

J *' j = J k (k 2 + m 2 Y ((q - k) 2 + m 2 y ' (R15) 

In addition to (F.6) we also need: 

(-2TdU°) 2 ab = - J2 U ab (U ae + U be ) + J2 U acU cb = ~ J2 & a b (K e + B' be ) + Y2Y. & a,B' eb ■ 

e c e e 

Since (— 2TdU°) ab etc. contains only at least 2-replica sums, it can be dropped. We define: 

Nl b , cd = K dM + Nl M (F.16) 

Kdfic = \TJlAdhc + C(q)6 ad v b v c + v a v d (C(q) 2 (-2TdU°) bc + C {q)\-2T dU°f bc 

+T5 ad (Jl 2 (-2TdU°) bc + Jl 3 (-2TdU°) 2 bc ) + ^TJ q 2j2 (-2TdU°) ad (-2TdU°) bc 
+TJl 3 (-2TdU°) ad (-2TdU°)l) . (F.17) 

We obtain: 

N ldM = Sad&bcLlb + S adPa, bc + S bc P bM ~ U bc V a V d C{qf 

+\TJ{ 2 U ad U bc + (v a v d C(q) 3 - TJl 3 U ad ) Ypv u fc - U UU bf + U cf )] (FAS) 

f 

L ab = \TJl l+ TJl 2 J2u b f (F19) 

/ 

Kbc = C(q)v b v c - TJ q h2 U bc + TJl fl J2l U bfUc f - U bc (U bf + U cf )] - \TJ{ 2 U bc ^ U af (F20) 

/ / 

Kad = VaVdC(q) 2 J2 Ubf ~ \TJ{ 2 U ad J2 U b f ■ (F21) 

/ / 

This yields: 

2 ~ ~ 4 

N a b ,cd = $ad5 bc L q ab + 5 ad Q q abc + 5 bc Q q bM - -(B' bc v a v d + B' ad v b v c )C(q) 2 + -Jl 2 B' ad B' bc 

^(v a v d C(q) 3 - 2Jl 3 B' ad ) J2[B bf B' fc - B' bc (B' bf + B' cf )] 

f 

^- 2 (v b v c c(q? - 2JIA) J2fia f B' fd - BUKf + %)\ 

f 

= TJl + 2JI J2(B'af + B' bf ) (F.22) 

/ 

Ql, bc = v b v c (C(q) + ^C{qf E B' af ) - ^(TJ^ + 2 J{ 2 £ B' a 

f f 



+ 7 

4 



af) 
f f 
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+Pls Yl&Af ~ B' bc (B' bf + B'cf)} ■ (P-23) 

/ 

F.4 Final calculation of the matrix M 

We now perform the matrix multiplication 

M\ cd = J2 M abtef Ni f!Cd = -\[B" ah Ni h ^ - s ab J2 KKJ > ( R24 ) 

e/ 9 

where we have defined: 

■^ab,cd = ( N ab,cd + N ba,cd ~ N aa,cd ~ N bb,cd) = L ab {8 ad 8 bc + 5 bd 5 ac ) - Ka^acd ~ L^bcd 

+8 a dQ% c + 8bdQl ac - 5adQ 9 a , ac - 8 b dQ bfic + hcQ bM + 8 a cQl M ~ &acQ aM ~ 8 b cQ b>bd 

-(U bc v a v d + U ad v b v c + U ac v b v d + U bd v a v c - U ac v a v d - U ad v a v c - U bc v b v d - U bd v b v c )I% 
J rTJ 2 ^{U ad U bc + U bd U ac — U ad U a c — U bd U bc ) 

+(v a v d C(q) 3 - TJl 3 U ad ) Y^PbfUfc ~ U bc (U bf + U cf )\ 

f 
+ (v b v c C(q) 3 - TJ{JJ bc ) J^PafUfd - U ad (U af + U df )} 

f 
+(v b v d C(q) 3 - TJl 3 U bd ) Y} U -f U fc - U ac(Uaf + U cf )] 

f 

+(v a v c C(q) 3 - TJl 3 U ac ) Y} U bi U fd - U bd (U bf + U d} )\ 

f 

-{v a V d C{qf - TJl 3 U ad ) Y^PafUfc - U ac (U af + U c f)\ 

f 

-{v a v c C{qf - TJ{ 3 U ac ) Y\U*fUfd - U ad (U af + U d f)\ 

f 

-(v b v d C(q) 3 - TJ{ 3 U bd ) Y} U MUfc - U bc (U bf + U cf )] 

f 

-(v b v c C(q) 3 - TJl 3 U bc ) J^PbfUfd - U bd (U bf + U df )} . 

f 

Performing the matrix product yields the parameterization of the matrix M., where we have not yet fully 
expanded in sums, as (the g-dependence is implicit): 

2 1 

a ab = ~ 7Z 1 4o// r ( R25 ) 

1 l - T B ab L ab 

%a = Oi aa y D a jL aa 
f 
Dad = —CH aa B ad L ad + a aa 2_^ B a j-(Q atad — Qajd) 

f 
Zab OL ab IJ ab J-Jaa 
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tabd — OiabB ab (Q a! bd ~ Qa. 



ad I 



Tj// T 



,c,ad) 



v ac d — 0L aa \B ad (Q ddc — Qd,ac) + B ac (Q ccd — Q c 

- J2 KfTJUUadUfc + u fd Uac - u ad u ac - Uf d U fc ) 
f 

+ ^ B af( U fcV a Vd + U ad V f V c + U ac V f V d + U fd V a V c 
f 

- U ac v a v d - U ad v a v c - U fc v f v d - UfdVfV c )I% 

Wabc = 

9abcd = UabB'a b { ~ (U bc V a V d + U ad V b V c + U ac V b V d + U bd V a V c 

- U ac v a v d - U ad v a v c - U bc v b v d - U bd v b v c )I% 
+ TJ 2 ^{UadU bc + U bd U ac — U ad U ac — U bd U bc ) 
+ (v a v d II - TJ{ 3 U ad ) J2l U bfUfc - U bc (U bf + U cf )\ 

f 

+ (v b vJ! - TJl 3 U bc ) YflafUfd - U ad (U af + U df )\ 

f 
+ (v b v d II - TJl 3 U bd ) Y} U -f U fc ~ Uac(U af + U cf )\ 

f 

+ (v a vJl - TJl 3 U ac ) Y^PbfUfd - U UUb f + U df )\ 

f 

- (v a V d II - TJl 3 U ad ) YflafUfc ~ U ac (U af + U cf )\ 

f 

- (v a V c I q 3 - Tj{ 3 U ac ) J^PafUfd - U UUaf + U d} )\ 

f 

- (v b v d I q 3 - TJ{ 3 U W ) Y^PbfUfc - U bc (U bf + U cf )} 

f 

- (v b vJi - TJl 3 U bc ) YpbfUf d - U bd (U bf + U#)]} . 

f 

We now finish the expansion in sums, using (F.9), (F.l 1) and defining the notations: 



A q 

A ab 



Al 



1 - *B ab Jl 
2 



l-AB"{2TI 1 )Jl 1 
*ab = -hK h + 2{Al h f{hl 2 B" ab £(#«, + & bf ) + nX b E °<*f)] 



L ab = TJl 1 + 2Jl 2 J2(B' af + B' bf ) 
f 



f 



f 



(F.26) 
(F.27) 
(F28) 
(F29) 
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G 




O = B (x) 
• = B(v2) 







G 




Figure 6: Some typical l/7V-diagrams and the classes they belong to. The wiggly line indicates the derivative, the black circle 
is a B (or derivative), the grey circle a B (or derivative). We have restrained from drawing a contribution due to the explicit 
f-dependence. These terms do not play a role in the argument, and are only tedious supplements one has to keep track of. 



h The result is given in the text in section 5, as well as: 

,4 

T 



P lVad = 2Jl 1 {A q aa)\^JhB" aa £ Kf + JlXa £ O aaf )^ ad + ^A^Xd £(K/ + &«) 



f 



C{q) 



+J{iA q aa B'" d J2 °adf ~ T (l-AB"(2TI 1 )J q 11 ) ^ B " a ^ Va ~~ Vf ^ d 



f 
2 ,„ ■=,, 



p lZab = ^iM^)% J iX £(#«/ + & u) + J lA £ °^s)Kt 



+7j;A ab Jl 2 b' ab 2_^ B' af + J\ x A q ah B ah 2_^ Oabf ■ 



(F.30) 



(F.31) 



G More remarks on the graphical method 

G.1 Diagrammatics 

The diagrammatic 1/JV expansion can be constructed by using the 
Theorem: 



B'(v 2 



d d 



dx dv 2 



B( X ) + ^BW(x(v),v 2 ) + -±- 2 B^\x(v),v 2 ) + 



(G.l) 
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Figure 7: Diagrams at order 1/N (first), 1/N 2 (next two) and 1/iV 3 (rest). 

£(!) ( x ) = ^ all 1PI diagrams with 1 loop (G.2) 

B {2) (x) = Y a11 1PI dia g rams wim 2 loops (G.3) 

... = ... (G.4) 

X (v) = v 2 + 2Th + 4I 2 [d x B( X (v 2 ),v 2 ) - d x B( X (0), 0)] . (G.5) 

Some explication and precisions are in order: 1 -particle-irreducible diagrams (1PI) are w.r.t. lines being 
correlators (vv), and vertices being B^(x)- The r.h.s. of Eq. (G.l) are diagrams made out of bare vertices. 
We have separated the ^-dependence from the explicit v -dependence: The latter are v's which are connected 
with a line. These are the terms in our 1/iV-calculation, which explicitly contain v. Note that v 's always 
pair. Side-chains only come from the fact that finally one inserts x- Note that x as defined here is an object 
which contains terms at all orders in 1/N. The diagrams are 1PI, a fact which is important for the order 
1/N 2 . It means that B^ does not contain the diagram made out of 2 closed loops, connected by a single 
line. 
Proof: 

Draw a collection of diagrams contributing to B (see figure 6). This drawing contains vertices made 
out of derivatives of B{v 2 ) (not B(x) - we have drawn the completely expanded diagram). Now derive 
that object with respect to v 2 , giving a couple of terms. Any of these terms singles out one B, namely the 
one derived. 

This B may be part of a tree, by which we mean that either it is a point or by cutting off one of the 
attached legs, the diagram will fall apart. Then it is contained in the first term on the r.h.s. of Eq. (G.l), 
since any attachment which can be made to it in the form of a tree, is taking care of by choosing the above 
given x- Note that for this to be true, x has to be exactly the object given above, i.e. on the r.h.s. of Eq. 
(G.5) there has to be the full B to all orders in 1/N. In the diagrammatic language this is clear: Having a 
higher-order diagram and taking the derivative at one of the tree-like vertices, this diagram may still contain 
an arbitrary loop somewhere attached to the tree. 

The B which has been derived may as well be part of a closed loop. By this we mean that when we cut 
off all parts of the diagram which can be disconnect from our chosen one by one cut, there remains more 
than the vertex itself. This object is of higher connectivity; it can either be a loop (at order 1/N); it can be 
a diagram in the form of an 8 or a circle to which one has added an additional line between two arbitrary 
chosen points on it (at order 1/N 2 ). Higher order diagrams are given in figure 7. 

G.2 An alternative derivation of the /^-function (T = 0) 

We now give a general derivation of the /3-function to all orders in 1/N, following our results of the last 
section. This derivation is restricted to T = 0, since it is rather complicated at finite T. To make the 
derivation more transparent, and to avoid having to derive with respect to v 2 on the r.h.s. of Eq. (G.l), 
we introduce the auxiliary function B(v 2 ,u 2 ). By u 2 we shall denote a pair of background-fields that 
are non-trivially connected by a line of propagators, whereas v 2 denotes the background-field which is 
inserted into B, but which is not connected to any propagator. Note that this decomposition is unique; 
that the paring of w's is natural, and that deriving with respect to v 2 , but not u 2 can combinatorically be 
interpreted as choosing any vertex B and deriving it once. This object is thus better fitted for calculations 
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th<mB{v 2 ) = B{v 2 ,u 2 ) 
We now start from a 

d 



. However the latter object is of course the only one with a physical meaning. 



modified version of Eq. (G.l), namely 



d v 2 B(v 2 , u 2 ) = — (B( X ) + SB( X , u 2 )) 

5B(x,u 2 ) = 5B[B'(x(v)),B"(x(v)),...,u 2 } 



(G.6) 



^BV[B'(x(v)), B»(x(v)), ■■■} + ^ (2) [S'(XM), B"( X (v)), ...,u 2 } + ... (G.7) 



X = X(v) = X(v, u 2 ) =v 2 + 4/ 2 ( d v 2B(v 2 , u 2 ) - d v2 B(0, 0) 
1 



v 2 + 



d v 2b(v 2 ,u 2 )-d v 2b(0,0) 



(G.8) 



where we are a little bit sloppy with the notations, suppressing the argument u of x( v )- We define the 
dimensionless quantities 

b(v 2 ,u 2 )=4eI 2 B(v 2 ,u 2 ) = 4A d m- e B(v 2 ,u 2 ) 

b (x)=4A d m- e B( X ) 
5b(x,u 2 )=AA d m- e 5B(x,u 2 ). 

As in the main text, we use the notation i n (p) = -^^ = ^M_ The /^-function is 



(G.9) 
(G.10) 
(G.ll) 



-m 



_d_ 
dm 



d v 2b(v 2 ,u 2 ) = ed v 2b(v 2 ,u 2 ) - (4:A d m- e ) 



md d 
dm dx 



(B(x)+5B(x,u 2 )) 



ed v 2b(v 2 , u 2 ) + jj^ (6(x) + Sb( X , u 2 )) (-m ® X 



dm 



+6 



d_ 



A=l L 



ld_ 

Xdx 



Sb(Xb' , A6q, . . . ,u 2 



(G.12) 



Note that in the last equation, we have been a little bit sloppy with the notation. What this means is that 
having rescaled B to bo, the m-dependence of the integrals is canceled. Thus we can evaluate all integrals 
at m 2 = 1. The derivative w.r.t. A is easily understood as follows: Having a diagram with n + 1 vertices, 
the integrals scale like m~ ne . First, this accounts for the factor of e. Second, in order to get the right 
combinatorial factor of n instead of n + 1, one has to subtract one contribution, which is done by the factor 
of 1/ A in front of 5b. 

We need two more equations. First, starting from Eq. (G.8) and deriving w.r.t. m, we obtain (exact!) 



— 771- 



dx 

dm 



-m^-\ [d v 2b(v 2 ,u 2 )-d v 2b(0,0)]^± 



-m— ) d v 2b(v 2 ,u 2 ) -d v 2b(0,0) , (G.13) 
dm 



where an equality would suppose that due to dimensional reduction d m B'(0) = 0. 
Deriving Eq. (G.8) w.r.t. v 2 , we obtain (also exact) 

dx i 1 



dv 2 



-d 2 v2 b(v 2 y 



Deriving Eq. (G.6) by v 2 gives with the help of Eq. (G.14) 

d 2 V 2b(v 2 ,u 2 : 



— Hx)+Sb(x,u 2 )}^- 2 



dx 2 

-^ Hx) + Sb( X , u 2 )] (l + - e d 2 v2 b(v 2 , u 2 ) 



(G.14) 

(G.15) 
(G.16) 
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Therefrom we infer that (to all orders) 

d 2 



d 2 2 b(v 2 ,u 2 ^ 



dx 2 



1 + ±d 2 2 b(v 2 ,u 2 ) 



This equation can also be written as 



dx 



2 L/„.2 „.2\ 



Hx) + Sb( X ,u 2 )} dv2b(V ' U 



(G.17) 



(G.18) 



This procedure can be repeated to obtain 



d 3 
dx 3 

d X 4 



[b (x) + 5b( X ,u 2 )} 
[bo(x) + Sb( X ,u 2 )} 



d 3 2 b(v 2 ,u 2 ) 



(l + i}d 2 2 b(v 2 ,u 2 )) 3 

d A v2 b{v 2 ,u 2 ) + \ (d 4 v2 b(v 2 ,u 2 )d 2 v2 b(v 2 ,u 2 ) - 3 [d 3 2 b(v 2 ,u 2 )\ 



Eliminating ^ and #^[&o(x) + 5b(x, u 2 )} from Eq. (G.12), we obtain 



dm dx 2 



d 



-mj—d v 2b{y 2 i u 2 ) = ed v 2b(v 2 , u 2 ) + ., T , ; — - 

dm 1 + ±d 2 2 b(v 2 ,u 2 



(l + ld 2 2 b(v 2 ,u 2 )) 5 
G. 
d 2 2 b(v 2 ,u 2 ) I 



(G.19) 
(G.20) 



l(- m ^)(^b(v 2 y)-d v ib(0,0)) 



+6 



d_ 

dX 



A=l 



ld_ 

Xdx 



Sb(\b' , Xbl . . .) 



(G.21) 



We now take the limit of v 2 ,u 2 — > 0. We suppose that (— m-^) (d v 2b(v 2 ,u 2 ) — 9^6(0,0)) — ► in 



q2 fo{ v 2 u 2\ 

that limit. Further d 2 2 b(v 2 ,u 2 ) can either remain finite or diverge. However - — f„ 2 , ' ,. remains fi- 

v v ' ' c l+-a Jr,^ 

nite whatever d 2 2 b(v 2 ,u 2 ) will do. Supposing that this argument is indeed correct (are there additional 
IR-divergences?) the conclusion is that 

1 d 



-m—d v 2b(0, 0) = ed v 2b(0, 0) + lim e— 

dm u,v^o OX 



A=l 



Xdx 



5b(Xb , Xb , . . .) 



(G.22) 



The /5-function thus is equivalent to 

d 



-m^d v 2b(v 2 , u 2 ) = ed v 2b(v 2 , u 2 ) + \d v 2b(v 2 , u 2 ) - d v 2b(0, 0)]d 2 2 b(v 2 , u 2 ) 
dm 



+ (e + d 2 v2 b(v 2 y)) 



d_ 
dX 

d 
d 2 v2 b(v 2 ,u 2 ) lim — 



A=l 



ld_ 

Xdx 



Sb(Xb' ,Xb'o,...) 



A=l 



^5b(Xb ,Xb' ,...) 



(G.23) 



Using Eq. (G.14), this can also be written as 

d 



-m^d v 2b(v 2 , u 2 ) = ed v 2b(v 2 , u 2 ) + [d v 2b(v 2 , u 2 ) - dM, 0)}d 2 2 b(v 2 , u 2 ) 
dm 



d 

+e dx 



A=l 



d_ 

dv 2 



jSb(Xb' , Afe'o, . . .) 



d 
-d 2 v2 b(v 2 ,u 2 ) lim — 

v,u-*0 OA 



A=l 



ld_ 

Xdx 



Sb(Xb , Xb' , . . .) 



(G.24) 
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Integrating the latter equation over v 2 , we obtain 

-m^-b(v 2 , u 2 ) = eb(v 2 , u 2 ) + -d v 2b(v 2 , u 2 ) 2 - d v2 b(v 2 , u 2 )d v 2b(0, 0) + e^- 
dm 2 dX 



a=i L 



A 



6b(\b' ,\b'o, 



d 
-d v 2b(v 2 ,u 2 ) lim — 



A=l L 



1 d 

j—5b(\b ,\b ,...) 



(G.25) 



which of course has to be read at u 2 = v 2 . In a final step, we want to reintroduce proper quantities. Noting 
that 



d v 2b(v 2 ,u 2 ] 
(which does not need u), we obtain 



b'(v 2 )-d v 25b( X (v),v< 



(G.26) 



d 



l ,n 



m^-b(v 2 ) = eb(v 2 ) + -b'(y 2 ) 2 - b'(v 2 )b'(0) - [b\v 2 ) - b'(0)]d v 2Sb( X (v),v 2 ) 



dm 



+d v 25b(x(0), 0)b'(v 2 ) + - [d v 25b(x(v),v 2 )} - d v 25b(x(v),v 2 )d v 25b( X (0), 0) 



d 
+e dX 



A=l 



-r5b(\b' , Xb'o, . . .) 



d 



[b'(v 2 )-d v 2b( X (v),v 2 )] lim — 



«^o d\ 



A=l 



l^-5b(\b ,\b' ,...) 
Xdx 



(G.27) 



(Of course d v 2Sb(x(0), 0) means first to derive and then to put the arguments to 0.) Also note that 
d v 25b(x(0), 0) is not 0, at least at order T. The /3-function at order \/N therefore is 

d_ 
dm 



-m^-b{v 2 ) = eb(v 2 ) + \b\v 2 ) 2 - b'(v 2 )b'(0) 



1 / d 
+ NVdX 



A=l 



-bU(\b> ,\b'>,...) 



[b\v 2 )-b\0)]d v 2b (1 \x(v),v 2 



+b'(v 2 ){d u 2b^( X (0),u 2 ) - lim J- 

u=0 v,u->0 d\ 



A=l L 



ld_ 

Xdx 



b^(Xb ,Xb' ,...) 



+0{ m 



This might better be grouped as 



(G.28) 



d 



1 u l 



- m ^-b( v *) = e b(v 2 ) + -b'{v 2 ) 2 - b'(v 2 )b'(0) 



dm 



1 / d 
+ N{ e dX 



A=l 



\bW(Xb' ,Xb»,...) 



— b ( v ) lim 



d 



v,u^o dX 



A=l 



ld_ 

Xdx 



bU(Xb ,Xb' ,...) 



-[b'(v 2 )-b'(0)]d v 2b (1 \x(v),v 2 ) + b'(v 2 )d u 2b^(x(0),u 2 



u=0 



+ °'^ 



(G.29) 



A caveat is in order: The rescaling has to be done on the level of bare vertices, not on the level of renor- 
malized ones. That would give a wrong result. However the derivative w.r.t. v 2 can be taken in any 
formulation. 
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G.3 The case d — 

As one can see from our final result for the /^-function in (7.7), specified to d = 0, it is a polynom in b of 
finite order, since the denominators present in h x (p), see (6.13) are identical 1. Since this come as quite a 
surprise, we show here why this must be so; actually it is a quite general feature of the 1/iV-expansion of a 
renormalizable theory in d = 0. 

We start to warm up with the diagram 




(G.30) 



The leading and next-to-leading contributions (in b) are 







(G.31) 



Our strategy is to apply Bogoliubov's 7£.-operation, see e.g. [67-70] and to show that only the first two 
terms contribute. Three remarks are in order: 

1.) The 7£.-operation in the context of a 1/iV-expansion is maybe not entirely natural. However we have 
in the above diagrams the property, that the terms already encountered at \/N (and thus taken care of in 
the /^-function at leading order) are exactly the iterated 1-loop diagrams, thus the first order in 1/iV. 

2.) In order to extract the (3-f unction from 1Z applied to a diagram, we only have to derive (w.r.t. m) the 
diagrams in the boxes, since only those are counter- terms. 

3.) Applying 1Z to a functional of the bare b gives the result as a functional of b. Thus the contribution 



5 (3(b) to the /^-function is 5/3(b) 



-m 



dm 



7£[diagram(6 )]- Now 



*Q^G 




(G.32) 



which give as the contribution to (3 



60(b) 




(G.33) 



The second diagram gives 



K 







(G.34) 



In D — 0, the last two diagrams cancel, and 



W) 




(G.35) 
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This gives also the ratio 6 to be found in the explicit formula: The 3 from the first diagram has to be set 
in relation to the combinatorial factor of 2 for the second one and a factor of 1/e for the second, which 
together give a ratio of 6; finishing the test. 

We now proceed to higher orders: First we remark that for the chain with n members (here n — 3) one 
can show recursively that 



ncxyo] 



o 



o 



Thus having n diagrams in the chain and deriving the n-dependence gives 



md 

dm 



K [OCX>] = ns 



o 



o 



o 



-\ n— 1 



Of course, in d = the latter vanishes for n > 1. 

It is now easy to see that only finite order terms can appear, q.e.d. 



(G.36) 



(G.37) 



H Details of the calculation of the /^-function at finite T 

H.1 Integrals appearing in the T > correction to the effective action 

The following functions have been defined but not given in formula (6.17) in the main text. 

If 1 



g5{a x ) 

g&M 

97(a x ) 
#9(00) 

9io(a>x) 
9ii(a x ) 

9i2{a x ) 



1 1 f 



Ad j p 
1 1 
~2A~ d 



H(p)K(p) + 2*2(p)i3(p)/i*(p) 2 



A d 

1 1 

'2A~d 



8A d 



16A d 



16A d 



c{p)i2(p) 2 h(p)h x (pY 
h(p)h (p) 



i2{p)h{p)h Q (p) 



- ln(l - a x i 2 {pj) + i2(p)h x (p) + -i 2 {pfh x {pf 



i2{p) 2 h x (p) + i 2 {pfh x (pf 



MP?h x (p) + ^2{pfh x (pf 



7 = -4«i»2 



^is(ao) = -0^2 / i2(p) 2 h (p) 



(H.l) 

(H.2) 

(H.3) 

(H.4) 

(H.5) 

(H.6) 

(H.7) 

(H.8) 

(H.9) 

(H.10) 
(H.ll) 
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and we recall that a cutoff A/m is implicit in all these rescaled momentum integrals. 
These functions are not all independent. Indeed, defining 



lm,nM ■= 4" [ f(pMp) n h x ( P r (H.12) 

A d J p 

one easily establishes the recursion relations: 

m7m+l,n(fl) = a 2 d a lmn{a) 

7m+l,n+l(o) = -7m+l,n(o) - 7m,n(«) • (H.13) 

They allow to obtain these functions from derivatives of 

7(a) = ~ / ln (l " ai 2 (p))f(p) (H.14) 

Specializing first to f(p) = 1, one finds using (H.13) 



7i,i(a) = aj'(a) (H.15) 

71.2(a) = 7'(o)-i? (H.16) 

7i,3(a) = -(V(a) - i?) - 70.2(a) (H.17) 

a 

72.2(a) = a 2 7 " (a) (H.18) 

72.3(a) = a 7 "(a)- 7 '(a)+i 2 (H.19) 

72, 4 (a) = /(a) - -(V(a) - i 2 ) + 70.2(a) (H.20) 

a 

Thus the functions gio(a), On(a), 012(a) and 013(a) can all be obtained from derivatives of 7(a) with 
/(p) = 1. Similarly, 04(a) and 05(a) can be all obtained from derivatives of 7(a) with f(p) = c(p), and 
similarly for g 6 (a), g s (a) and g 9 (a) with f(p) = i 3 (p), g 7 (a) with f(p) = i 3 (p)c(p). 

One may further attempt to relate 7(a) for different functions f(p). Since d m 2l 2 = — 27 3 , one can use 
that 

*s(p) = ^h(p) - ^' 2 (p) - ^A9 A i 2 (p) (H.21) 

where A = A/m. Integration by part yields identities such as 

f % d ( p d p + Xd x )H(t 2 (p)) = -d [ % d H(i 2 (p)) + Xd, [ % d H(i 2 (p)) , (H.22) 

J V J V J V 

which can be used to relate the integrals. 

H.2 Calculation of the T > /^-function 

Below we compute —md m b(x) at fixed A/m, thus we truly compute —{md m + Ad\)b(x). It is therefore 
useful at T > only for d < 2 when all integrals are UV convergent and the limit A/m — > 00 can be 
taken with no further redefinitions. The calculation of the /3-function for d > 2 and T > requires further 
redefinitions and will eventually be detailed elsewhere. 
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Taking into account all m-dependence at T > in (6.17) one obtains (up to an additive constant): 

(H.23) 

- \-m.rP nJ + ^ \-mPP (ti(rr) - h'(0\\] + ^L(-m.ff a,] + ^L(-m.ft 

5a 7 



Pi[b](x) = -md m b(x) 

" L [-md° m a x ] + — 



8[V(x) - &'(0j] \-- md ^ h '^ ~ 6 '(°))] + ^ { - md ° ma0) + ^- md "mTm) 



+^(-«M - eh - (b> - 6'(0))(6' 1 - 61(0)) - T m i + ^ 

where b\ is given in (6.17). 
One uses that: 



T 

J. rr 



b'(x) b'((0) 
;i + 6"(0)/e) 2 e 



-md^a x = ea x + [b'(x) - b'(0)}a' x + 



T 

J. n 



1 + 6"(0)/e 



-m^(6'(x) - 1/(0)) = [e + 6"(x)][6'(x) - 6'(0)] + 

,2 



T 



— m<9° a = ae + T„ 



a 



l + 6"(0)/e 

(a — e) 2 a 



[6"(x) - 6"(0)] 



a -e 



(H.24) 
(H.25) 
(H.26) 



and —md m T m = —9T m . After a rather tedious calculation one obtains the form (7.13) given in the text 
with the following definitions (for d < 2): 



g 4 {a) 



-4#i(a) - e(e + % 4 (a) + e ag' 4 (a 



95(a) = -gi(a) - 2g 2 (a) - e6g 5 (a) + e ag 5 (a) 

g s (a) = -2g 2 (a) - e6g 6 (a) + e 2 ag' 6 (a) - eg 4 (a) 

2 
97(a) = -92(a) - 4513(a) - eg 5 (a) + e(e - 9)g 7 (a) + e 2 ag' 7 (a) 



98(a) 



ae 2 4 

— — [g 6 (a) + g$(a) + 0(^7(0) + 59(a))] + -51(a) + eg 4 (a) - e6g 8 (a) 



+a 



+a" 



— o9i(a) + -92(a) + eg 5 (a) + e(e - v)g 9 (a) + e g 8 (a 



,2 / , 



-,9i{a) 



2„/ , 



-g 2 (a) + -53(a) + e g g (a 



9io(a) 



(e + 20)g 1Q (a) - eg 4 (a) + e ag' 10 (a 



9n(a) = -2e 6gn(a) + 54(a) - 655(a) - eg 6 (a) + e ag' n (a) 



(H.27) 
(H.28) 
(H.29) 
(H.30) 

(H.31) 

(H.32) 

(H.33) 

(H.34) 
(H.35) 
(H.36) 



512(a) = e (e - 29)g 12 (a) + 55(a) - 657(a) + e ag 12 (a) 
0(a) = — — [5 6 (a) + 057(a) + g 8 (a) + eg 9 (a) + e 2 g' 10 (a) + ae 2 g' u (a) + a 2 e 2 g' 12 (a) + e 3 7 + ae 3 5i 3 (a)] 



+25 4 (a) + eg' 8 (a) - 2eV + a 
+,r e 3^ (o) _ MM + 2 ^ (o ) 



2 
e 2 (e - 2%i 3 (a) - -g 4 '(a) + 2g' 5 (a) + 2g' 6 (a) + 659(a) 



ij,(a) = e[g" (a) + ag'i^a) + a 2 g" 2 (a) + eai 3 (a) + aeg' 13 (a)] 
$(a) = (e - a) 2 [5' 10 (a) + a^ u (o) + a 2 ffi 2 (a) + e 7 + eo^is(o)] 



(H.37) 

(H.38) 
(H.39) 
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I Integrals 



1.1 Definitions 



h = A d 



1 



k (k 2 + m 2 ) n 



A d : 



2r(3 - d/2) 

(47f) d / 2 



The momentum dependent ones are 

hip) ■- 

Up) ■- 
hip) ■-- 

Dimensionless rescaled variants: 



k (k + p/2) 2 + m 2 ik- p/2) 2 + m 2 

1 1 

,, [(k + p/2) 2 + m 2 } 2 (k - p/2) 2 + m 2 

1 1 

k \{k + p/2) 2 + m 2 ] 2 [(k - p/2) 2 + m 2 ] 2 ' 



«n(p) : = 



Up) 



A, 



m=l 



b n .- 



A, 



171=1 



1.2 Integrals in fixed dimensions, general formulas 

The general case can be treated as follows: 



Jn,m\P) ■- 



1 



AiJk [ik-l) 2 + l] m [ik+l) 2 + l\ n 
Using the usual Schwinger-parameter representation, this can be written as 

1 ( 1 



Jn,m\P) 



a/3 „2 



a 



r(n)r(m) \A d Jk j Ja0>o 

Tin + m-d/2) f ^ m ~ x 

2r(3-d/2)r(n)r(m) Jfro 



1 /r- 1 (a + /3)"3 e - (a+/3) e"^ y 



;i+/3) 



n+m 



1 + T^W P . 



We make the change of variables (5 = ^£- 

Tin + m-d/2) r ' 



Jn,m\P) 



21 2~ m ~ n 



. dss m - 1 il-s) n - 1 \l + sil-s)p 2 ] 
2r(3 - d/2)r(n)r(m) 7 

r(n + m-d/2) 

T(3 - rf/2)r(n)r(m)2«+ m 

- 1 dy (l+y^^) m - 1 (l-y^^) w - 1 + (l+ v ^^)"- 1 (l-v^^) m - 1 

o V^y 2 



(1.1) 

(1.2) 
(1.3) 

(1.4) 
(1.5) 
(1.6) 

(1.7) 
(1.8) 



(1.9) 



(1.10) 



1+2/ 



p- 



(I.H) 
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where we have used another new variable y = 4s (1 — s). Note that the large fraction ^ in the above 
expression is such that only integer powers of (1 — y) survive. Some simplifications occur for n — m, and 

hip) =hM- 



JnM 



T(2n - d/2) 



r(3-d/2)r(n) 2 2^J ^T=y- 



dy -.y n - 1 



1+y 



P 



2T|-2n 



hip) = JiM = o 



dy 



^=V 






(1.12) 
(1.13) 



1.3 d = 

h = hi?) = 1/4 , A d = 4 



e = 4 



(1.14) 



1.4 d = 1 



A 

H 

hip) 

hip) 

hip) 



4 

2 



3(4 + p 2 ) 

12 +p 2 

3(4 +p 2 ) 2 

2(20 +p 2 ) 



3(4 + p : 



2^3 • 



(1.15) 
(1.16) 
(1.17) 

(1.18) 

(1.19) 



arctanh 



hip) 



2+p 2 



IpIv^ + p 2 



1.5 d = 2 

1 



A, 



2^ 



(1.20) 



In (2+p 2 + \p\ ^/4 + p 2 ) - In (2 + p 2 - |p| ^4+"^) 2 arcsinh( 



IpI- 



2|p|v^T7 



«3(P) 



+ 



2arcsinh(§) 



8 + 2 P |p|(4 + p 2 



U(p) = « 2^1(2,3,5/2,-^/4). 




bl\/4 + p 2 

(1.21) 

(1.22) 



(1.23) 



1.6 d = 3 

1 



A, 



i 2 (p) = 2 



arctan(^- 



IpI 



8tt 



[ln(2-i|p|)-ln(2 + i|p|)] 



(1.24) 

(1.25) 
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hip) 

u(p) = 



p 2 + 4 



(p 2 + 4) 2 • 



(1.26) 

(1.27) 



57 



J Summary of Notation 



symbol 


definition 


defined in equation 


e 


e = 4-d 




CO 


( = roughness, 9 = d — 2 + 2( (thermal exponent) 




u(x), v(x) 


uix) = field, v(x) = u(x)/vN 




In 


r ._ / i r _ a m ' e a _2r(3-rf/2) 

~ Jki^ + w?)^ z a e ' a (4vr)^ 


(2.14) 


Up) 


I (v) - f l l 




(4.18) 


^ yt,> ' J k ik+p/2) 2 + m 2 (k-p/2) 2 + m 2 


hip) 


Iiv)- f l l 


> 


(4.20) 


' Jk [(k + p/2) 2 + m 2 ] 2 (k - p/2) 2 + m< 


hip) 


Iiv)- f l l 


(4.21) 


' Jk [(k + p/2) 2 + m 2 } 2 [(k - p/2) 2 + m 2 } 2 


C(p) 


C(p) : =(p 2 + m 2 )- 1 




l"n\P)i ^n 


in( P ) := J f } 

Ad 


h 

' %n '~ A., 


m=l 


(6.15), (6.16) 


cip) 


c{p):={l+p 2 Y l 




B(...) 


second cumulant of bare disorder 




B(...) 


second cumulant of renormalized disorder (not rescaled) 




B' ab , B" ah , B' ab , etc. 


B' ab := B'(xab), B'i h := B"(xab), Kb ■= B\v 2 ab ), etc. 




b(z) 


b(z) := AA d m^- e B(zm- 2 ^) 


(6.10) 


Xabix), Xabix) 


auxiliary fields 




Xabix) 


Xabix) := Xabix) + Xbaix) ~ Xaaix) - Xbbix) 




~ab ^-.ab 

X.V » X.V 


Xv = Xab\X) \v(x)=v> Xv = Xv >^\~n) 




Xv Xv Xab 


xf:=v 2 ab + 2Th + 4I 2 


B 'ab ~ \ \B'aa + B bb j 






H v ip) 


H h) ■ B "^ 


) 


(4.27) 


v[Ph l-4hip)B»ixv 


HM 


H (v) ■ & ' {X) 




(6.8) 


^ v " 7 ' l + A[I 2 -I 2 ip)]B"ix) 


Kip) 


h H . b%x) 




(6.13) 


x{P) ' l + [i 2 -i 2 ip)]b»ix) 


T 


temperature 




T 


T m := ATA d m e /e 
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